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ABSTRACT 

Using //5T/COS/STIS and HIRES/Keck high-resolution spectra, we have studied a remarkable Hi ab- 
sorbing complex at z = 0.672 toward the quasar Q1317+277. The Hi absorption has a velocity spread of 
Av = 1600 km s" 1 , comprises 21 Voigt profile components, and resides at an impact parameter of D = 58 kpc 
from a bright, high mass (logM v i r /M Q ~ 13.7) elliptical galaxy that is deduced to have a 6 Gyr old, solar 
metallicity stellar population. Ionization models suggest the majority of the structure is cold gas surrounding a 
shock heated cloud that is kinematically adjacent to a multi-phase group of clouds with detected C III, Civ and 
Ovi absorption, suggestive of a conductive interface near the shock. The deduced metallicities are consistent 
with the moderate in situ enrichment relative to the levels observed in the z ~ 3 Lya forest. We interpret the 
Hi complex as a metal -poor filamentary structure being shock heated as it accretes into the halo of the galaxy. 
The data support the scenario of an early formation period (z > 4) in which the galaxy was presumably fed by 
cold-mode gas accretion that was later quenched via virial shocking by the hot halo such that, by intermediate 
redshift, the cold filamentary accreting gas is continuing to be disrupted by shock heating. Thus, continued fil- 
amentary accretion is being mixed into the hot halo, indicating that the star formation of the galaxy will likely 
remain quenched. To date, the galaxy and the Hi absorption complex provide some of the most compelling 
observational data supporting the theoretical picture in which accretion is virial shocked in the hot coronal 
halos of high mass galaxies. 
Subject headings: quasars: absorption lines 
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1. INTRODUCTION 

It is well accepted that galaxies are intimately linked to 
the gaseous cosmic web and that the evolution of galax- 
ies is governed in large part by the dissipative response 
of baryonic gas due to the trade off of cooling and dy- 
namical timescales as it accretes into dark matter ha- 
los (e.g.. IBinnevI \Wf% i Rees & Ostrikefl \Wf% Islffl [19771 
White & Rees 1978). In general, distinct modes of accretion 
are believed to operate and the mode is dependent primar- 
ily up on the dark matter halo mass (e.g., jBirnboim & Dekell 
2003; Dekel & Bimboim 2006), with some dependence on 
environment (e.g., iKeres et al.l 120051) and on feedback (e.g., 
Ivan de Voort et alj|201 llT ~ 

"Hot-mode" accretion is the mechanism in which inflowing 
gas is shock heated as it is compressed by the hot hydrostatic 
gas halo. This mode dominates around high mass galaxies, 
where the dynamical time is shorter than the cooling time. 
For the most part, the gas accretes into the halo, but not nec- 
essarily onto the galaxy itself. "Cold-mode" accretion is the 
mechanism primarily around low mass galaxies, where the 
cooling time is shorter than the dynamical time so that hot hy- 
drostatic halos do not form and cold gas accretes directly onto 
the galaxy. At lower redshifts, as densities decrease, the cool- 
ing time is generally longer than dynamical time and the rate 
of infalling material decreases so that cold mode accretion is 
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a minor channel of accretion o nto galaxies (e.g.. IKeres et alj 
12001 iDekel & Birnboimll2006T) . though it remains an impor- 
tant channel for the gro wth of galaxies at all redshifts (e.g., 
Ivan de Voort et alj|201 ll) . 

A possible third mode is the case in which cold dense 
filaments can penetrate directly into a hot halo of a mas- 
sive galaxy. If the conditions are favorable for short cool- 
ing times in the filament, then shock heating can be avoided 
and t he filamentary gas can directly accrete onto the galax y 
(e.g JPekel & Birnboiml2006tlvan de Voort & Schavel201 ll) . 
However, though cold streams can penetrate the hot at- 
mospheres of massive halos at z > 2, th is process signifi- 
cantly diminishes at lower r edshift (e.g.. IKeres et al.l 120091 : 
iFaucher-Giguere et al.l 120 lTb . though there is plausible ev- 
idence of this proce ss occuring in some z < 1 galaxies 
dKacprzak et alj|201 UlRibaudo et alj201 ltlThom et alj|201 It 
iKacprzak et alj|2012l) . 

One observable signature of a filament may be a "complex" 
of Hi absorption with a large velocity spread (Av > 1000 
kms" 1 ). Alternatively, such Hi absorbing complexes may 
trace the warm hot intergalactic medium (WHIM), or the in- 
tracluster and/or intragroup medium. If an observed H I ab- 
sorption complex arises from a filament, it is plausible that 
the filament may be accreting from the intergalactic medium 
into a galaxy halo or even directly onto a galaxy. Thus, Hi ab- 
sorption complexes provide unique astrophysical laboratories 
for placing constraints on our understanding of the intergalac- 
tic medium, and the processes giving rise to extended galaxy 
halos in the context of dark matter overdensities. 

At high redshifts (z > 2), kinem atically extended H i ab- 
sorbing complexes were studied by ICowie et~aT1 (11996b with 
Keck/HIRES spectra. They deduced that these structures 
were filamentary in nature. At lower redshift {z < 0.22) 
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Shull et alJ (fl998h.lTripp et al.1 d200lh . lShull et al.1 d2003l) . and 
Aracil et all (l2006al) studied three different Hi complexes 
with high resolution ultraviolet spectra in which metal lines 
were detected. For these three complexes, several galaxies 
were found in the vicinity, suggesting moderate size groups. 
Each of these studies favored a different scenario of expla- 
nation for the physical picture of the H I complex, including 
O VI arising in "nearside/backsi de" shocked infall into the po- 
tential well of the galaxy group (S hull et ail l2b03). intragroup 
gas or an unyiria lized filamentary structure through the group 
(Tripp et al.l l2.001l) . and tidally stripped material from one of 
the nearby galaxies (lAracil et al.ll20 06a). 

The spectrum of the quasar Ql 3 17+277 (TON 153, CSO 
0873, J131956+272808, V = 16.0, z em = 1.017) exhibits a dra- 
matic Av > 1000 km s" 1 Hi complex observed in Lya ab- 
sorption at z = 0.672. Several optical and ultraviolet spec- 
tra of the quasar are available and have been the focus of 
this Hi co mplex and/or the Lyman limit metal-line system at 
z = 0.660 dSteidel & Sargenllll992tlBahcall et al.lfT993l | 1996t 
Churchill et al.ll2000l: iDing et al.ll2005t IChurchill etaill2007l: 
Kacpr zak et alJ l201lt IChurchill et alJ 12012; Ka cprzak et alJ 
20121) . Besides being at intermediate redshift, what is unique 
to this Hi complex is that it lies at 58 kpc p rojected from a 
single bright elliptical galaxy at z ga i = 0.6719 (IChurchill et al.1 
120071) . 

An analysis of the z = 0.672 Hi complex was presented 
in IChurchill et al.l (120071 hereafter Paper I) based upon HST 
G160L/G190HFOS spectra (PID 2424, PI: J. N. Bahcall), an 
HST E230M/STIS spectrum (PID 8672, PI: Churchill), a nd a 
Keck/HIRES spectrum of the quasar (see lChurchillll 19971) . In 
the FOS spectrum, the Hi complex was found to comprise 
five components of optically thin Hi absorbing gas, which 
span a velocity range of Av = 1400 km s" 1 . No metal-lines 
were clearly detected in the FOS, STIS, and HIRES spectra. 

Galaxies at the absorber re dshifts were first reported by 
ISteidel. Dickinson. & Perssonl (119941) as pail of the their Mg II 
absorption-selected galaxy survey. The quantified morpho- 
logical and spectral properties of these galaxies were pre- 
sented in Paper I based upon an HSTPNFPC2 F702W image 
(PID 5984; PI: Steidel) of the quasar field and Keck/LRIS 
spectra o f the galaxies. Updated analysis of the galaxies is in- 
cluded inlKacprzak et aTT(l201 ll) . IChurchill etall (120121) . and 
iKacprzak et alJ (120121) . In summary, the galaxy at the redshift 
of the Hi complex has z ga i = 0.6719 and an impact parameter 
of D = 58 kpc. It is classified as a late-type E/S0 galaxy. The 
galaxy at the redshift of the z = 0.660 absorber, which is a rich 
metal-line Lyman limit system, has z g - d \ = 0.6610, an impact 
parameter of D = 104, and is an inclined Sab galaxy. 

In Paper I, we deduced that the gas in the Hi complex has 
column densities in the range 14.5 < logA^(Hl) < 15.5, tem- 
peratures in the range 5.0 < logT < 5.5, and upper limits on 
metallicity in the range -1.0 < logZ/Z Q < -3.0. We fur- 
ther deduced the complex is consistent with a combination of 
photo and collisi onal ionized gas. Based upon expectations of 
simulations (e.g. lDave et alll999h . we favored a shock heated 
structure with chemical enrichment consistent with the high 
redshift intergalactic medium. Our interpretation was that the 
Hi complex is a shock heated filamentary structure that orig- 
inated as a photoionized diffuse phase of gas at high redshift 
and is accreting in the vicinity of the z s - d \ = 0.6719 elliptical 
galaxy. 

However, the FOS data did not provide the information 
required to examine the detailed structure of the kinematic 



and ionization conditions. Thus, we were unable to con- 
strain the kinematic relationships between the putative shock 
heated gas, the photo and/or collisionally ionized gas, and 
the Zgai = 0.6719 galaxy. Such information is critical for 
examining the physics of gas accretion onto galaxies and 
for comparing with other observa tions and w i th cos molog- 
ical simulations. For example, ICowie etail (fl996) claim 
that clustered Lya lines at z = 3 are observed to have 
higher ionization conditions at the velocity extremes and 
suggest this layered structure is a signature of collapsing 
structures. While this may be a signature for intergalac- 
tic filaments, it is not clear that accretion of a filament into 
a galaxy potential will yield the same velocity-ionization 
structure. In theoretical t reatm e nt (e.g.. | Birnboim & Dekell 
120031: iDekel & Birnboiml 120061; iBirnboim et al 1 120071) and 
cosmological simulations ("e.g.. |Keres et al.l |2005l 120091 : 
Ivan de Voort et al.ll201 lllvan de Voort & Schayal201 ll) accre- 
tion onto massive galaxies is expected to shock heat upon 
entry into a hot coronal halo, and the resulting kinematic - 
ionization structure may reflect this very different process. 

In order to better study the Q13 17+277 Hi complex at 
z = 0.672, we have obtained a high resolution R = 18,000 
HST/COS spectrum of the quasar, with focus on the H I Ly- 
man series lines and Ovi AA1031, 1037 absorption. We fur- 
ther improved our knowledge of the star formation history, 
age, and mass of the z ga i = 0.6719 galaxy by obtaining multi- 
band imaging of the quasar field. Our motivations include 
(1) a thorough examination of the kinematic, ionization, and 
chemical structure of the H I complex, (2) a direct comparison 
with the galaxy properties to help place the galaxy-absorber 
system in the context of galaxy evolution scenarios predicted 
by theories and cosmological simulations for an improved 
interpretation, and (3) examination whether there is a con- 
nection between the z = 0.672 Hi complex and the galaxy- 
absorber pair at z = 0.661, perhaps in the form of "bridge" of 
weak Ly a absorption between the galaxies. 

The paper is structured as follows: Reduction and analysis 
of the imaging and spectroscopic data are presented in § [2] 
Our analysis of the galaxy images is presented in §|3]and our 
analysis of the absorption line data is presented in § @] In § [5] 
and §0 we present our ionization modeling and resulting con- 
straints on the physical conditions of the absorbing gas. We 
discuss and interpret the data and our findings in § [7] We pro- 
vide brief a conclusion in § [8] Throughout this work, we as- 
sume the cosmological parameters fl\ = 0.7, VL m = 0.3, VL^ = 0, 
and ft = # /(100 km s" 1 Mpc" 1 ) = 0.7 (based on lJarosik et all 
I20T1 . 

2. OBSERVATIONS, REDUCTIONS, AND CALIBRATIONS 

We have acquired new data on the z = 0.660 and z = 0.672 
absorption and the two galaxies Gl and G2. First, we ob- 
tained a cycle-17 R ~ 18,000 COS spectrum (PID 11667; 
PI: Churchill) of the quasar covering the transitions examined 
in the FOS spectra. Our goals include measuring higher de- 
tection sensitivities and detailed kinematics of the absorption 
lines. Second, we obtained g', r', /', and K s ground-based im- 
ages of the quasar field. The multi-band images provide col- 
ors from which galaxy stellar populations, metallicities, and 
masses can be estimated. 

2.1. Ground-based Images 

As part of a larger campaign for imaging Mg II absorption 
galaxies, we obtained g', r', i' band images of the quasar field 
using the Seaver Prototype Imaging camera (SPIcam) on the 



Quenched Cold Accretion 



3 



(a) g' band 



G2 



(b) r' band 



* Gl 



G2 



(c) i' band 



Gl 



G2 



(d) Ks band 



• Gl 



G2 • 



(e) 






• Gl 


• / 




*- 






G2 








(f) 



* 

Gl D=58kpc 




6250 6500 6750 7000 7250 7500 7750 8000 



(g) 



G2D=104kpc 




5750 6000 6250 6500 6750 7000 7250 7500 7750 



FIG. 1 . — Images of the quasar field and spectra of the absorbing galaxies. — (a-c) The 43" X 60" sections of the ground-based APO/SPIcam images centered 
on the quasar through the SDSS g', r' , i' niters. Galaxies Gl and G2 are identified. — (d) A 43" X 60" section of the KPNO/IRAM NICMOS image through 
the K s filter. — (e) A 12" X 22" section of the //S77WFPC2 image. Galaxies Gl and G2 are labeled. Also identified and labeled is object G3, a compact 
galaxy (with unknown redshift) at angular separation 2.1" in the north by north-east direction from the quasar. — (f) A 4" X 4" view of galaxy Gl at D = 58 
kpc from the quasar. — (g) A 4" X 4" view of galaxy G2 at D = 104 kpc from the quasar. — (h) The Keck/LRIS spectrum (counts versus wavelength) of Gl, 
yielded z ga ] = 0.6719 from the Call H&K doublet. — (i) The Keck/LRIS spectrum (counts versus wavelength) of G2, yielded z gii i = 0.6610 from the [Oil] A3727 
emission line. 



Apache Point Observatory's (APO) 3.5-m telescope. The de- 
tector is a 2048 x 2048 CCD with 24 fim pixels, giving an 
unbinned plate scale of 0.14" pixel" 1 and a field of view of 
4.78' x 4.78'. We binned the CCD 2x2 during readout for a 
plate scale of 0.28" pixel" 1 . Multiple g', r 1 , and i' band im- 
ages were obtained on the nights of 2006 March 24-25, for 
which the seeing varied between 1.2-1.5", and additional g 1 
band images were obtained the night of 2007 March 15 for 
which the seeing was 0.7". The total summed exposure times 
are 5190, 4630, and 4350 seconds for the g', r', and i' filters, 
respectively. 

Each frame was reduced using standard IDL and IRAF0 
scripts and tasks. Flat fielding incorporated a combination 
of dome and sky flats. Cosmic rays were removed in each 
individual frame. The astrometry was calibrated by position 
matching of USNO A2.0 stars in the field. Final images were 
obtained by coadding the individual calibrated frames. 

The photometric zero points were obtained using stars in 
SDSS images. Color terms are required because the filter 
plus detector throughput of the APO facility is not identical 

6 IRAF is written and supported by the IRAF programming group at 
the National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NOAO is operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc. under cooperative agreement with the National Science 
Foundation. 



to that of the SDSS facility. These color terms, which are of 
order 0.1, were determined from SPIcam and SDSS images 
of roughly 30 quasar fields from our more extensive database. 

The near-infrared band images were obtained on 1994 
February 24 using the Kitt Peak Mayall 4-m telescope through 
the K s filter (1.99-2.32 fim) with the IRIM NICMOS III 
256 x 256 array camera. These images were obtained as 
part of the campaign culminating in the work of Steid el et al.l 
( 1994). The field of view is 154" x 154" with plate scale of 
0.6" pixel" 1 . The NICMOS images were reduced using the 
contributed IRAF package DIMSUM3 Updated photomet- 
ric zero points were determined using star s from the 2MASS 
point-source catalog (Strutskie et al. 2006|). 

All p hotometric measuremen ts were conducted using SEx- 
tractor (iBertin & Arnoutsl 1996b: we adopted the AUTOMAG 
results. The dust maps of SchlegiEitir|C998) were used to 
correct for Galactic dust extinction. 

In Figures QJ-QJ/, we present 43" x 60" sections of the 
ground-based images centered on the quasar Q1317+277. A 
12" x 22" section of the //S77WFPC2 F702W image is pre- 
sented as Figure QJ. Galaxies Gl and G2 labeled in all im- 
ages. In the HST image, note the object to the north by 



7 http://irafMoao.edu/imf/flp/contrib/ctimsumV3/. DIMSUM was con- 
tributed by P. Eisenhardt, M. Dickinson, S. A. Stanford, & F. Valdez. 
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north-east within 2.1" of the quasar. We have no estimate 
of the redshift of this object, which we label G3. Reduction 
and analysis of the F7 02 W image was described in Paper I, 
iKacprzak etaD (120111) . and lChurchill et all (|2012l) . Expanded 
views of the two galaxies are presented in Figures [T]f andQJ. 
The Keck/LRIS spectra of galaxies Gl and G2 (originally de- 
scribed in Paper I) are presented in Figures \T}i andQj. The Gl 
galaxy redshift was determined using Gaussian centroiding to 
the Call absorption features and the redshift of G2 was deter- 
mined using Gaussian centroiding to the [Oil] A3727 emis- 
sion line. 

2.2. HST/COS spectrum 

We obtained a cycle- 17 R ~ 18,000 COS spectrum (PID 
11667; PI: Churchill) of the quasar covering the transitions 
first examined in the FOS spectra. Two NUV/G185M spec- 
tra were obtained 2010 May 26 and optimally co-added. The 
first was centered at 1921 A for a 5420 sec exposure, and 
the second was centered at 1941 A for a 4970 sec exposure. 
The overlap region was 2223-2037 A on Stripe C, which pro- 
vided a total of 10,390 sec of integration on the the Hi ab- 
sorption complex. The FUV/G160M spectrum was obtained 
2010 June 26 centered at 1600 A for an 12,580 sec exposure. 

The s pectra were reduce d following the procedures out- 
lined by IShaw et al.l (120091) . We continuum fit the lower or- 
der shape of the spectrum using the IRAF sfit task and refined 
the higher order contin uum features using our own code Fitter 
(ChurchilfetaDHQOi. 

For the z = 0.660 and z = 0.672 absorbers, the Ly /3 and Ly 7 
absorption lines were captured in the NUV on Segment A, the 
Ly 8 was not captured (fell between segments), and the higher 
order Lyman series lines were captured on Segment B. 

3. IMAGE ANALYSIS: GALAXY PROPERTIES 

A re-analysis0 of the //S77WFPC2 F702W image was un- 
dertak en, presented, and fully described by Churchill et al . 
d2012l) . We adopt the measured quantities from that work. 

The Gl quasar-galaxy impact parameter is D = 58 kpc. 
The galaxy photometric properties are Mb = -21.6, Mk = 
-23.0, and B-K = 1.4 (AB). The B-band luminosity is 
Lb /Lg = 1.28, where we use Mg from the fit with redshift re- 
ported by iFaber et"aT1 (12007b . From an analysis using GIM2D 
( Sima rd et al.l2002l) . we measure a half-light radius of r/, = 4.2 
kpc, disk scale length of = 0.9 kpc, and bulge-to-total ratio 
of B/T = 0.99. The galaxy classifies as an E/S0 based upon 
its C-A morphology (Abraha m et al.l[l996l) . The galaxy in- 
clination is i = 15.9°, and the angle between the quasar sight 
line and the major-axis of the projected ellipse of the galaxy is 
= 22. 1 °. For the properties of ga laxy G2. see | Churchill et all 
(2012|)and our companion paper (IKacprzak et alj 120121) . 

Photometric analysis of the ground-based images yielded 
dust-, color-, and seeing-corrected AB apparent magnitudes 
for galaxy Gl of m g > = 23.4 ± 0.03, m r > = 22.2 ± 0.03 m r = 
21.01 ±0.03 and m^ s = 19.4±0.1. A color composite image 
of the ground-based images can be viewed in our companion 
paper (Kacprza k~et al.l2012l) . from which it can be ascertained 
that galaxy Gl is clearly redder than galaxy G2 and that the 
other galaxies in field (see Figures QJ-Q}/) are likely at sub- 
stantially different redshifts. 

For galaxy G3, we measured an F702W apparent magni- 
tude of m = 23 .4 ±0.2. This value is based upon a dithered 

8 We note that an incorrect ^-correction resulted in an overestimate of Mb 
and Lb jL* B in Paper I. We also converted all magnitudes to the AB system. 



TABLE 1 
Galaxy Properties 



Property 


Gl 


G2 a 


G3 b 


M»/M 


3 x 10" 


1 x 10" 


1 x 10 10 


M VU /M Q 


5 x 10 13 


8 x 10 12 


8 x 10" 


Tylr, K 


1 x 10 7 


3 x 10 6 


7 x 10 5 


v circ > km s 


550 


280 


150 


Rvir, kpc 


750 


380 


180 


M gas /M Q 


2 x 10 10 


1 x 10 10 


7 x 10 9 


M baIy /Mo 


3.2 x 10" 


1.1 x 10" 


1.7 x 10 10 


ft 


0.05 


0.09 


0.34 


" Taken from Kacprzak et a 


1.(2013). 





b Assuming z = 0.672. 



co-added image constructed by A. Shapley, in which she re- 
moved the quasar via point spread function (PSF) subtraction. 
The quoted uncertainty is statistical based upon the sky back- 
ground; due to the PSF subtraction, the error could be sub- 
stantially underestimated. Assuming galaxy G3 is at the red- 
shift of the Hi absorption complex, i.e., z — 0.672, and adopt- 
ing the quoted apparent magnitude, we compute Mb = -19.3 
(L B /L B = 0.16) and M r = -20.7 (SDSS r-band). At this red- 
shift, the impact parameter to G3 is D = 14.6 kpc. 

In Figure [2] we plot the g'—r 1 colo r s vers u s i'—K s colors 
for Gl and G2. Following iBell et al.l ( 120031) . iFontana et ail 
(2004), and lSwindle et al.l (1201 11) . we used stellar population 
models to determine stellar masses, M*, of galaxies Gl and 
G2 from the o bserved colors. We employ ed the stella r popula- 
tion m odels o flBruzual & Char iot (2003) assuming a Chabrier 
(2003) initial mass function and an exponential star forma- 
tion history with an e-folding time of 1 Gyfl Also shown in 
Figure [2] ar e the locust of observed c olors as a function red- 
shift for the lBruzual & Charlotl (120031) stellar population mod- 
els for the metallicities [Z/H] = -0.4, 0.0, and +0.4. We find 
galaxy Gl is consistent with a ~ 5.8 Gyr old, solar metallicity 
stellar population with a formation epoch of z = 4. The stellar 
population models also yield the galaxy Zf-band mass to light 
ratio, (M/L) K . From this ratio and the /^-band magnitude, we 
estimate logM*/M© = log(M / ' L) K -0A(M K -M K . Q ) = 11.5 for 
galaxy Gl, where Mk.q = -3.28 is the solar value. 

Using the technique of halo abundance match- 
ing, the galaxy virial mass, M^ u , can be es ti mated 
from the stellar mass (e.g. , IConrov & Wechslerl 1 2009 ; 
[ Behroozi. Conrov. & Wechslerl 120101 Ovloster et al.l I2OIO ; 
IStewardfeOl ll) . Abundance matching assumes a monotonic 
functional relation between M v ; r and by assigning the 
number of halos with Mhaio > Myit equal to the number 
of galaxies with M„. d \ > M*. As such, it matches the halo 
mass and stellar mass functions globally with a roughly 0.25 
dex uncertainty in M v ; r at fixed M t , primarily due to the 
systematics in estimates of ( Behro ozi et al.l l20ip|). We 
emplo yed the parameterized functions presented by (IStewartl 

mob. 

For galaxy Gl, we obtained M vi l = lO 13 - 7 " 14 - 4 M fi 
where the l ower value is given by the IConrov & Wechslerl 
(2009) and iMoster et all (120101) fi ts and the higher value 
is from the iBehroozi et all (120101) fi{3 We adopted the 

9 The stellar population models were generated using the web service EZ- 
GAL (www.baryons.org/ezgal). 

10 Abundance matching is the most accurate for M v j r = 10 11-13 Mq. For 
less massive halos, the stellar mass function is not tightly constrained. For 
higher mass halos, particularly massive ellipticals, there is substantial scatter 
between the published abundance matching predictions/relations. Thus, for a 
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IConrov & Wechslerl (120091) and iMoster et al.l (120101) values 
M vir = 10 13 7 M n . 




g-r 

FIG. 2. — The g'—r' versus i'—K s color-color diagram showing the measured 
dust-c orrected colors of galaxy Gl and G2. The curves are Bruzual & Chariot 
120031) stellar population models for metallicities [Z/H] = -0.4, 0.0, and +0.4 
assu ming an exponen tial star formation rate with an e-folding time of 1 Gyr, 
and a|Chabrier (2003) initial mass function. The arrow provides the redden- 
ing vector for E(B—V) = 0.1 in the rest frame of the galaxies. 

Under the assumption that galaxy G3 resides at z = 
0.672, we use the ACDM Bolshoi Simulation Database of 
Truiill cTGomez et ail d201 ll) to estimate a virial mass of M vir = 
10 M Q using abundance matching. This mass is the aver- 
age of 18,500 galaxies in the absolute magnitude bin —21 .7 < 
M r < -19.7 with average M r = -20.7. The abundance match- 
ing in this database is constrained by the observed luminosity- 
circular velocity relation, the baryonic Tully-Fisher relation, 
and the circular velocity function, allowing all types of galax- 
ies to be included. Using the par ameterized abunda nce match- 
ing function for M*(M vil ,z) of ( IMoster et al.ll2010t) . we esti- 
mate that G3 has a stellar mass of M* ~ 10 10 2 M Q . 

The average gas mass, M gas , of a galaxy with stellar mass 
M* can be es timated using the parameterized relation of 
iStewarj (1201 l j) based upon th e baryonic Tully-Fisher rela- 
tion study of McGaugh (2005 b, the stellar, gas, and dynam- 
ical mass relation of Erb et al. (2006), and galaxy gas fraction 
stellar mass study of IStewart et al.l (120091) . We obtain M gas = 
10 io.2±o.2 Mq for galaxy G1 rn] and estimate Mgas = io' ^j 

M© for G3. Thus, we deduce averaged baryonic gas fractions, 
f g = M gas /(M* +M gas ), of 5% for Gl and 34% for G3. 
The virial radii, virial temperatures, and circular velocities 

single case, abundance matching is not highly robust for mapping M v ; r from 
M* when M v j r > 10 13 Mq. (Stewart 2012, private communication) 

" Estimates of the gas mass is based upon the M gas -M* correlation for 
disk galaxies, and is therefore not directly applicable to elliptical galaxies. 
However, as shown in Fig. 2 of[Stewart et al. (2009), forM, > 10 11 Mq, the 
gas-poor disk galaxy gas frac tions reasonably match those of the red galaxy 
sample of Kannappan 1 2004). Thus, we can crudely apply the relation to 
M* > 10" Mq elliptical galaxies ( Stewart 2012, private communication). 



for galaxies Gl and G3 w ere computed using the relations of 
iBrvan & Normanl (11998b for fi R = 0. For galaxy Gl, we ob- 
tained /? vir = 750 kpc, r vil = 1 x 10 7 K, and v cilc = 550 km s -1 
and for galaxy G3, we obtained R v i r = 180 kpc, T v [ r = 7 x 10 5 
K, and v c ; lc = 150 km s" 1 . In Table Q] we summarize the de- 
duced galaxy properties, which are representative averages for 
galaxies with their observed photometric properties. 

4. SPECTRAL ANALYSIS: ABSORPTION PROPERTIES 

In the FOS spectrum, a very broad H I absorption complex 
was observed (first reported by Bahcall et all 11993b . In Pa- 
per I we fitted the Ly a with five Gaussian components at red- 
shifts z = 0.66914, 0.67157, 0.67355, 0.67559, and 0.67707, 
respectively. The total rest-frame equivalent width was de- 
termined to be Wf(Lya) = 2.87 A over a rest-frame veloc- 
ity spread of 1420 kms" 1 . From the STIS spectrum, the 
3 a Civ A 1548 equivalent width limit was estimated to be 
W r (1548) < 0.03 A at z = 0.6719. In the FOS spectrum, the 
3 a Ovi A 1031 equivalent width limit was estimated to be 
W r (1031) < 0.21 A. For Mgll A2796 in the HIRES spectrum, 
we obtained W, (2796) < 7 mA to 3 a. These values assumed 
unresolved absorption. 

For this work, we objectively locate absorption lines or 
place limits on th eir equivalent widths empl oying the opti- 
mized m ethods of Schnei der et al.l (119931) and lChurchill et ail 
(1999b) as modified using the m ethods for unresolve d lines 
and pattern noise developed by lLawton et al.l (12008). We 
adopt a 5 a detection threshold and quote 3 a limits. 

4.1. Neutral Hydrogen Lines 

In Figure [3] we present the Lya, Ly/3, Ly7, Lye, and Ly £ 
absorption observed in the COS spectrum as a function of 
rest-frame velocity relative to the z gai = 0.6719 galaxy. No Hi 
absorption was detected for Lyman series lines higher than 
Ly£. There is significant blending in the Ly/3, Ly7, Lye, and 
Ly £ lines, which are identified in Figure|3] when possible. 

The curves through the data are \ 2 minimized Voigt pro- 
file (VP) fits obtained using our code MINFIT (IChurchilll 
fl99"7t IChurchill & Vogtl I200U IChurchill et al.l 120031) . The 
ticks above the normalized continuum provide the VP com- 
ponent velocities. During the fitting, the COS instrumental 
line spread function (ISF) was convolved with the VP model. 
The COS ISF appropriate for the spectrograph settings and 
the observed wavelength of each transition was determined 
via interpolation of the on-line tabulated data (cf jDixon et al.l 
2010; lKrissll20TTh . 

For the VP fitting, pixels compromised by blending were 
masked out of the least-squares fit vector. We present the 
results of our VP modeling in Table [2] We fitted a total of 
21 components, or "clouds" (MINFIT returns the minimum 
number of components based upon their statistical signifi- 
cance through a series of F-tests and confidence level checks). 
Assuming thermal broadening, we converted the Doppler b 
parameter into the "cloud" temperature. 

The strong Lya absorption centered at v ~ -530 kms" 1 
(z — 0.6690 ) is blended with Ly/3 a t z = 0.9776 associa ted 
with the Ly a absorption identified bv lBahcall et al.l ([1996) at 
z = 0.9778 in the FOS spectrum, which we confirmed in the 
STIS spectrum. Prior to VP fitting the H I complex, we de- 
blended the v ~ -530 km s" 1 Ly a feature using VP modeling 
and employing both the COS and STIS ISFs. In Figure gjz, 
we present the VP fits at z = 0.97763. The flux decrement 
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FIG. 3. — Voigt profile decomposition of the Hi complex, including Ly/3, 
Ly7, Ly£, and Lye (Ly<5 was not captured). The vertical dashed line is the 
redshift of galaxy Gl. The component parameters centroids are marked by 
ticks above the continuum and the fitted component parameters are listed in 
Tables[2]and[3] Several blends are present and identified as: (1,2) O VI A1031 
associated with the z = 0.6610 metal system; (3,4) Cm A977 at z = 0.6610; 
(5) Ly <5 at z = 0.7738; (6) C III A977 from clouds #7 & #8 of the H I complex; 
(7) Ly^ z = 0.5349; (8) Galactic Civ A1550; (9) LyC at z = 0.6610. Those 
marked with "(?)" may well be Lya lines in that there is no corroborating 
data to suggest they are metal lines nor higher order Lyman series lines. 

of the Ly/3 VP model was then subtracted from that of the 
v ~ -530 km s" 1 Lya absorption from the Hi complex. The 
resulting deblended profile is illustrated in Figure SI?. 

Prior to and following the deblending process, we experi- 
enced difficulty obtaining a satisfactory simultaneous fit to 
the Lya and Ly/3 absorption in the v ~ -530 km s" 1 feature. 
Though the deblending provided an improved VP model, we 
caution that the number of VP components and the result- 
ing column densities and Doppler b parameters we adopted 
should be viewed with discretion (see comments for Table[2). 

4.2. Metal Lines 

We examined the COS, STIS, and HIRES quasar spec- 
tra for associated metal-line transitions, including the 
Mgn AA2796,2803, Siiv AA1393, 1402, Civ AA1548, 1550, 



TABLE 2 

Voigt Profile Decomposition of the Hi Complex 



Cld# 




V 1 :i 
v ct 


loeMHlI 


b(H I) 


T 






[Km s j 




[Km s j 


Liu jy| 


l b 


0.668645 


-582.1 


13.24 ±0.59 


31.9±31.3 


62.0 ±60.9 


2 b 


U.UUOOU 1 


—543 5 




14 7 ± 9 9 


13 2 ± 8 9 


ii> 
j 


u.ooyuuu 


— 3 lo.O 


1/1 ^1 1 f\ 11 


n o 1 i n "7 
IZ.o ± 1U. / 


1U.U =t o.4 


4' 


0.669112 


-498.5 


14.51 ±0.14 


3.2±2.1 


0.6 ±0.4 


5 b 


0.669199 


-482.9 


14.07 ±0.20 


4.0±1.8 


1.0 ±0.4 


o 




JZU.O 


1 1 78 4-0^8 

Ij.jO nz u.jo 


1 ^ S -I- 1 1 


117-1- in s 

1*+. / TX 1U. J 
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0.671203 


-123.8 


14.13±0.10 


11.5 ±2.9 


8.1±2.0 


8 


0.671244 


-116.6 


13.69±0.14 


43.3 ±14.4 


114.5 ± 38.1 


9 


0.671617 


-49.6 


12.93 ±0.94 


5.8=b 15.5 


2.1 ±5.5 


10 


0.671739 


-27.7 


13.70 ±0.20 


18.3±7.7 


20.5 ±8.6 


11 


0.672319 


+76.1 


13.43 ±0.19 


88.6 ±49.2 


479.3 ±266.1 


12 


0.672994 


+197.0 


13.59 ±0.28 


26.4 ±11.9 


42.5 ±19.2 


13 


0.673291 


+250.3 


13.93 ±0.14 


27.6 ±10.9 


46.4 ±18.3 


14 


0.673562 


+289.9 


13.38 ±0.26 


10.5 ±6.9 


6.7 ±4.4 


15 


0.673787 


+339.1 


13.68=1=0.15 


38.4±14.1 


89.9 ±33.0 


16 


0.674297 


+430.5 


13.11 =1=0.12 


26.1 ±10.1 


41.6± 16.1 


17 


0.674872 


+533.5 


13.70±0.13 


34.6 ±8.8 


73.0±18.6 


18 


0.675272 


+605.1 


14.26 ±0.05 


33.7 ±4.2 


69.4 ±8.6 


19 


0.675789 


+698.8 


13.82 ±0.04 


29.0 ±3.4 


51.5±5.9 


20 


0.676366 


+801.2 


13.10±0.12 


20.5 ±7.9 


25.7 ±9.9 


21 


0.676776 


+874.7 


13.86 ±0.04 


33.5 ±3.6 


68.3 ±7.3 



a Velocities are measured with respect to z„ a i = 0.6719. 

b See text for discussion of these components, which required a deblending treat- 
ment. We also fitted this region with a single VP component at v = —516.9, with 
log TV = 14.88 ± 0.02, and b = 32.76 ± 0.861 km s~' . However, the VP model 
provided a poor match the structure in the core of the Ly a feature. 



z =h S = 0.9776302 z IbI = 0.6733566 
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FIG. 4.— (a) The Ly/3 line (COS) from the z a i, s = 0.9778 Lya absorber 
(center panel) is blended with the v ~ —530 km s~' absorption of the Hi 
complex. The solid curves are a two-component Voigt profile model of the 
Hi absorption centered at z a b* = 0.97763, including the constraint from Lya 
(STIS) and Ly7 (COS). — (b) Subtraction of the z abs = 0.9776 Ly/3 Voigt 
profile model from the affected component of the H I absorption complex. 

and Ovi AA103 1,1037 zero-volt resonance doublets, and 
transitions from Sill, Sim, Cll, Cm, etc. In Figure |5] 
we present the wavelength regions of the quasar spectra 
corresponding to detected metal lines Cm A977 (COS), 
Civ AA1548, 1550 (STIS), andOvi AA1031, 1037 (COS). We 
also show the Ly a absorption over the same velocity range. 
VP components #7 through #1 1 are shown, where the "cloud" 
number is given in column 1 of Table [2] Only clouds #7, #8, 
and #10 have detected metals. 

In Paper I, we reported no detected metal absorption in 
Civ (STIS) and Ovi (FOS). However, we clearly detected 
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FIG. 5. — Expanded region of the Voigt profile decomposition of the de- 
tected metal lines over the redshift window for clouds #7 through #11 of the 
H I complex. The vertical dashed line is the redshift of galaxy G 1 . 

Ovi absorption in the COS spectrum in the velocity range 
-200 < V < kms" 1 relative to z gal = 0.6719. We also 
detected Cm A977 absoiption at v ~ -120 kms" 1 . We 
then re-examined the C IV absorption, and formally detected 
Civ A 1548 aligned in velocity with the Ovi and Cm absorp- 
tion. The Civ A1550 absorption, which we had interpreted as 
noise in Paper I, is detected at the 3.5 a level. 

In Figure [5] the curves through the data are fitted VP com- 
ponents. For these fits, we modified MINFIT to hold the 
VP component redshifts (velocities) and Doppler b param- 



eters constant, allowing only the column densities to be \ 2 
minimized, and fitted the metal lines using the velocities 
and Doppler b parameters from the VP fits to the Hi Ly- 
man series lines. As with the Hi fits de scribed above, the 
COS ISF (iDixon et alJ l20Tot iKrissI [20TTh and the STIS ISF 
(iBostroem et al.l l2010h appropriate for the spectrograph set- 
tings and the observed wavelength of each transition was de- 
termined via interpolation of the on-line tabulated data. 

When a VP component column density was objectively de- 
termined to be insignificant, MINFIT returned an upper limit. 
Several tests for significance were conducted during the least- 
square fitting convergence; in short, for a column density 
to be deemed significant, the fitted VP value and its uncer- 
tainty had to be inconsistent with the upper limit on the col- 
umn density (in that precise region of the spectrum and for 
that Doppl er b parameter), and th e integrated apparent optical 
depth ( Savage & Sembach 1991) had to be consistent with the 
VP component and inconsistent with the the upper limit on the 
column density. 

In Table [3] we list the metal-line absorption properties 
for all clouds. In MINFIT, the equivalent width limits are 
determined directly fr om the spectra using the m ethods of 
ISchneider etai] (119931) and lChurchill et al.l (Il999bl). but mod- 
ified f or partially or fully resolved features (see Lawton et al.1 
2008). For each transition, we use the Doppler b parameter 
from the VP fit to the H I series convolved with the appropri- 
ate ISF to measure the upper limit on the equivalent width in 
each cloud. The uncertainty in the Doppler b width is used 
to determine the 1 a spread in this limit. The column den- 
sity limits and the 1 a spread in these are determined from the 
curve of growth. 

5. IONIZATION MODELING 

We employed our own photo+collisional ionization code 
dChurchill & Klimekll2012h. whi ch is very similar to the code 
LINESPEC dVerner & Iakovlevll 19901) . The code treats pho- 
toionization, Auger ionization, direct collisional ionization, 
excitation-autoionization, photo-recombination, high and low 
temperature dielectronic recombination, charge transfer ion- 
ization by H + , and charge transfer recombination by H° and 
He . Metals up to zinc can be incorporated, and all ioniza- 
tion stages for each elemental species are modeled. The code 
is appropriate for optically thin gas in which no ionization 
structure is present. 

The input cloud parameters are the hydrogen number den- 
sity, « H , equilibrium kinetic temperature, T, and the mass frac- 
tion of metals, Z. Solar abundan ce mass frac tions are taken 
from T able 1.4 of Draine] d20Tl . based upon lAsplund et al.1 
(2009), though the code has the option o f modi fying the 
relative abundances. A lHaardt & Madaul (1201 ll) ionizing 
spectrum is used for the ultraviolet background (UVB). 
However, a stellar ionizing spectrum or a combined stel- 
lar plus UVB ionizing sp ectrum can be incorporated (see 
iChurchill & Klimekll2012l for details). 

The code obtains an initial guess solution for the density of 
each ionic species based upon the assumption of adjacent ion 
stage ionization and recombination balance (i.e., neglecting 
Auger and charge transfer processes). Then, the rate matrix is 
solved using the code dqed.f, a Hanson/Krogh nonlinear least 
squares algorithm w ith linear cons traints based on quadratic- 
tensor local model (lHansonlll986l) . The outputs of the ion- 
ization code are the electron density, and the ionization and 
recombination rate coefficients, ionization fractions and the 
number densities for all ionic species. 
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TABLE 3 

Selected Metal-Line Measurements 
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a Velocities are measured with respect to Zgai = 0.6719. 

Since the galaxy Gl is 58 kpc from the location where 
the gas is probed in absorption, and the stellar population is 
clearly dominated by red stars, we assume a UVB-only ioniz- 
ing spectrum. We also assume a solar abundance pattern and 
include metals up to iron (however, we omitted lithium, beryl- 
lium, boron, fluorine, and the nobel gas elements for which 
many of the ionization and recombination rates are not well 
determined). 

For a given cloud model, the resulting column density of 
ionic species X is obtained by 



log /V(X) = log /V(H i) - log/,,, + log n x - log n„ , 



(1) 



where 7Y(H I) is the H I column density from the VP fit to the 
data, «„ is the input hydrogen density of the cloud model, 
and where f m is the ionization fraction of H° and n x is the 
number density of ionic species X output by the code, re- 
spectively. Note that the metallicity of the model is implied. 
In practice, n x in a given model scales directly in propor- 
tion to Z. Thus, if the metallicity of the cloud model is 
Z m , for which the resulting density of species X is n x m , then 
log« x = log« x ,m+(logZ-logZ m ). For a given n H , such scaling 
of the models is valid only if /„, is independent of metallicity, 
which holds for log Z/Z Q < — 1. 

Comparing Eq.[T|to the measured values of 7V(X) from our 
VP fitting, we can constrain the metallicity and n„ of each 
cloud. The range in these constrained quantities are based 
upon the 1 a uncertainties in the measured column densi- 
ties. Only three clouds have detected metal lines and there- 
fore measured N{X). These clouds are #7, #8, and #10. We 
examine these clouds in the following subsections. 

5.1. Clouds #7 and #8 



Clouds #7 and #8 are located at v ~ -120 km s 1 with re- 
spect to the systemic velocity of Gl and are separated by 
Av = 7 km s" 1 , which is on the order of half a single resolu- 
tion element. The Hi profile is highly suggestive of a narrow 
plus a broad component combining to yield a deep core with 
a broad wing. The VP components reflect this profile shape, 
with cloud #7 providing the profile core, log 7V(Hl) = 14.1 and 
£>(Hl) = 8 km s" 1 , and with cloud #8 providing the broadened 
wings, logTY(Hl) = 13.7 and b(Hl) = 43 km s" 1 . The require- 
ment for the broader cloud is significant at the 99% confidence 
level, even though some contribution to the blue portion of 
the broad wing is due to the extremely broad VP component 
(cloud #1 1) at v = +76 km s" 1 . 

At v ~ -120 km s" 1 , Cm, Civ and Ovi are formally de- 
tected in the spectra. However, from the VP fitting, both 
clouds #7 and #8 have associated Cm and Ovi, whereas only 
cloud #7 has associated Civ. Inspection of Figure [5] shows 
that the C IV absorption appears to be offset in velocity from 
the Ovi. However, the Civ is measured in the STIS spec- 
trum, whereas the Ovi and Cm are measured in the COS 
spectrum. Thus, the measured velocity offset relies heavily 
on the accuracy of the wavelength zero points between the 
two independent spectra. 

Assuming thermal broadening, the implied temperature of 
cloud #7 is 8000 K, which is too cool to host strong Ovi 
and Civ absorption. The implied temperature of cloud #8 
is T= 1 14, 000 ±38, 000 K, which is consistent with the tem- 
perature at which C IV absorption peaks for collisional ion- 
ization equilibrium. To analyze this absorption structure and 
constrain its hydrogen density and metallicity, we opted to 
assume a sing le ph ase of gas that is dominated by cloud #8 
(however, in § !6.5l we also model these clouds as two distinct 
ionization phases). We added the Hi and metal column den- 
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FIG. 6. — (a) Ionization models of clouds #7+8 for a metal mass fraction of logZ/Zo = —2.0 with solar abundance ratios and temperature T = 1 14, 000 ±38, 000 
K. The solid curves are the model column densities of Ovi (red), Civ (green), and Cm (blue) and the dotted curves accounting for the uncertainties in the 
temperature. The measured values of N(0 VI), N(C IV) and N(C III) are provided in the panel to the left. — (b) The logZ-log n H plane showing the constraints on 
the models from the measured data (same color scheme). Solid black curves provide the overlap regions. 
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FIG. 7. — (a) Ionization model of cloud #10 for a metal mass fraction of logZ/Zg) = —1.4 with solar abundance ratios and temperature T = 20, 000 K. The solid 
curves and data points are as described in Figure|6] whereas the dotted curves account for a 0.2 dex error in logiV(Hl). — (b) The constrained range of the metal 
mass fraction on the logZ-logfiH plane, accounting for the column density limits on the carbon ions. 

tion is from O vi, with logyV(Ovi) = 13.6 ±0.1. 

A high abundant the +5 ion in a T = 20,000 K cloud may 
be somewhat surprising. For n H < -3 and log Z/Z Q =-2, anal- 
ysis of the rate coefficients output by the ionization code show 
that the balance of O 44 , O" 1 " 5 , and O" 1 " 6 is dominated by pho- 
toionization and recombination via charge exchange with H + , 
which dominates over free electron recombination by 2 orders 
of magnitude for O" 1 " 5 . 

For cloud #10, the measured upper hmit on Cm is 
logyV(Cin) < 12.7 and on Civ islogyV(Civ) < 13.2. As with 
with our exercise for clouds #7+8, we found that the column 
density limits for Mg II, Si ill, and Si IV did not contribute to 
constraining the gas properties. 

In Figure|7] we illustrate our analysis of cloud #10. The ion- 
ization code was run for -7 < log« H <— 1 for T = 20,000 K. In 
Figure [TJz, we show the column densities for Cm, Civ, and 
OVI for logyV(Hl) = 13.7 ±0.2 and metallicity logZ/Z Q = 
-1.4. The metallicity and hydrogen density constraints for 
cloud #10 are shown in Figure 13?. 

If we account for the uncertainty in the implied temperature, 
we find that for T = 28,000 K, the curves in Figure [7J? are 
unchanged for logn H > -5.3, and are shifted upward in the 
diagram by no more than 0.7 dex (higher Z at a given n H ) for 
logn H < -5.3. As temperature is decreased, the curves shift 



sities from clouds #7 and #8. For this exercise, we found that 
the column density limits for Mg II, Si III, and Si IV did not 
contribute to constraining the gas properties. 

In Figure [6] we illustrate our analysis of clouds #7+8. 
The ionization code was run for -7 < log«„ < — 1 for T = 
114,000 K, and for T = 152,000 and T = 76,000, which 
brackets the 1 a uncertainty in the temperature. In Figure^, 
we illustrate the model column densities for Cm, Civ, and 
Ovi for logyV(Hl) = 14.3 and metallicity logZ/Z Q = -2.0. In 
Figure [6f?, we show the metallicity constraints as a function 
of « H for each metal ion, where the range in the constraints 
are based on the uncertainty in the temperature and the uncer- 
tainty in the Hi and metal ion column densities. The region of 
overlap indicates where the three ions provide consistent con- 
straints on « H and log Z/Z Q . Allowing for the uncertainties in 
the data, the gas is constrained to have -4.2 < logn H < -3.7 
and -2.0 < logZ/Z Q < -1.6 by all three metal ions. Note 
that metallicity is contrained primarily by 7V(C III), whereas 
hydrogen density is constrained primarily by N(Ovi). 

5.2. Cloud #10 

Cloud #10 is located at v ~ -30 km s" 1 with respect to the 
systemic velocity of Gl. It is optically thin in neutral hydro- 
gen, with logTV(Hl) = 13.7 ± 0.2 and an implied temperature 
of T = 20, 500 ± 8600 K. The only detected metal line absorp- 
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to the upper right in the diagram such that at T = 12. 000 K, a 
minimum metallicity in the allowed range -2.3 < log Z/Zq < 
—1.8 occurs at logn„ > -3.7 instead of -2.5 < logZ/Z Q < 
-2.0 at logn H = -5.3. It seems less likely that Ovi absorption 
would be as strong as we detect in 12,000 degree gas. 

5.3. Clouds with No Detected Metal Lines 

The remainder of the H I absorption complex is character- 
ized by multiple components with gas temperatures in the 
range 20,000-90,000 K, though a hot, nearly million degree 
component is also present (cloud #11). For the clouds with 
no detected metals, the densities, ionization corrections, and 
therefore metallicities c anno t be constrained without addi- 
tional assumptions. In § 16.51 we assume hydrodynamic equi- 
librium and undertake an analysis of the clouds with only up- 
per limits on the metal line column densities. Below, we mo- 
tivate the assumption of hydrodynamic equilibrium as a rea- 
sonable scenario. 

6. INFERRED PHYSICAL CONDITIONS OF THE GAS 
6.1. Ionization Equilibrium and Cloud Stability 

In general, ionization equilibrium is valid when the colli- 
sional ionization time scale, t C0 u, and the photoionization time 
scale, T p h, are shorter than the cooling time of the gas, r CO oi- 

For a monatomic gas, the cooling time is the ratio of the en- 
ergy per unit volume and the energy loss per unit volume per 
unit time, t coo i = Q/Q = 3 / 2 (n e + n A ) £7 /n e it A A(T,Z), where n e 
is the electron density, n A is the total density of all ions, and 
A(T,Z) is the specific cooling function for metallicity Z. We 
obtain n e and n A from our ionization code as a function of n H 
and T. Our cloud models are optically thin, such that n A oc n„ 
and n e oc « H at a given T, we thus find that the behavior of the 
cooling time follows r coo i oc n~ l at fixed T. 

We adopt the Z = specific cooling function of 
ISutherland & Dopitd tfT!^) for a gas in collisional ionization 
equilibrium (also see Dopita & Sutherland 2003). Since our 
models include phot oionization, we corrected fo r photoion- 
ization heating (see Osterbrock & Ferland 2006), which re- 
duces the magnitude of Q, and therefore increases the esti- 
mated cooling time. We found that photoionization heating 
from the UVB at z = 0.67 is negligible for the density and 
temperatures ranges we explored. 

Metals strongly contribute to the cooling rate, so the Z = 
curve provides an upper limit on the estimated cooling time 
under the assumption of ionization equilibrium. However, 
the effect is negligible for log Z/Z Q = -3 and is a maximum 
difference of 0.4 dex at l og T ~ 5.3 for logZ/Z Q = -2 (see 
ISutherland & Dopitalll993l) . 

For a given ion, the collisional time scale is well approxi- 
mated as Toon — 1 /[n e (a KC + a co \{)], whereas the photoioniza- 
tion time scale is r p h ~ 1 /R p h(Jv), where a rec (T) and a co n(r) 
are the recombination and collisional ionization rate coeffi- 
cients, and R p h{J^) is the photoionization rate for the given 
ion, respectively. We obtain these quantities directly from our 
ionization code. Since n e oc « H at a given T, we have t co h oc 
at fixed T. Thus, the ratio t coo \/t co \\ is effectively constant as 
a function of n H for fixed temperature. Note that for lower 
temperatures, a co \\(T) vanishes so that r co n becomes equal to 
the recombination time scale, r rec — 1 /n e a lec . 

6.2. Equilibrium Based On Hydrogen 

In what follows, we briefly explore selected clouds using 
representative n H densities for the purpose of illustration. In 



Table[4] we present a summary of cloud models for the n H val- 
ues log« H = -5, -4, and -3, where ion specific quantities ap- 
ply for hydrogen. The hydrogen photoionization rate from the 
UVB at z = 0.67 is R ph = 3.0 x 10" 12 s _1 , yielding r ph = 1 x 10 4 
yr. As we shall see, the photoionization time is significantly 
shorter than the cooling times for all models; as such, only the 
longer collisional time scales constrain whether the clouds are 
in ionization equilibrium. 

For clouds #7+8, at T = 1 14,000 K, a rec = 6.2 x 10" 14 and 
atoll = 5.3 x 10~ 9 cm 3 s" 1 , respectively. Assuming logn H = -4, 
we obtain r coo i ~ 3 x 10 s yr and r co u = 5 x 10 4 yr. With 
Tcooi/Vcoii — 5 x 10 3 , hydrogen is clearly in ionization equilib- 
rium in clouds #7+8; the short photoionization and collisional 
ionization time scales relative to the recombination time scale 
are indicative of the highly ionized condition (f m = 2 x 10~ 6 ). 
The condition r p h < r co n indicates photoionization marginally 
dominates over collisional ionization. 

For cloud #10, at T = 20,000 K, a rec = 2.4 x 10" 13 and 
atoll = 4.3 x 10~ 12 cm 3 s" 1 . This cloud is constrained to have 
logn H < -4; assuming log«„ = -5, we obtain r coo i ~ 1 x 10 8 
yr and t co ii = 6 x 10 8 yr. For this assumed density, we find 
t coo i/t co u ~ 0.2, suggesting that ionization equilibrium may 
be marginal in cloud #10 for hydrogen (the gas tempera- 
ture evolves on a similar time scale that the balance can be 
achieved). As listed in Table |4] the thermal time scale is 
shorter than the ionization time scale for the presented range 
of n H . 

For cloud #11, the hot, log T = 5 .68 K, cloud, we have a rec = 
1.7 x 10" 14 and a coU = 2.7 x 10" 8 cm 3 s" 1 . For logn H = -4 
and -5, we obtain t coo i — 7 x 10 9 " 10 yr, and t co i] = 1 x 10 4 " 5 
yr, where the longer times correspond to the lower density. 
Cloud #11 has t coo \/t co u ~ 7 x 10 5 , indicating an ionization 
equilibrium condition. Note that if logn„ = -4, we have r p h ~ 
T co n, indicating that hydrogen ionization (/ HI = 3 x 10~ 7 ) is 
driven by both photo and collisional processes. If log n H = -5, 
then / HI = 6 x 10~ 8 , which is dominated by photoionization 
(Tcoii/Tph = 10); the larger ratio due to the longer recombina- 
tion time scale results in a higher ionization condition for this 
lower density. 

With regard to the remaining clouds in the Hi complex 
with only limits on the metal line measurements, we selected 
clouds #15 [logJV(Hl)= 13.7,logr=4.95],#18 [log/V(Hl) = 
14 . 3 , log T = 4 . 84] , and #20 [log 7Y(H I) = 1 3 . 1 , log T = 4 .4 1 ] as 
representative. Of clouds with no detected metals, #15 has the 
highest temperature and an intermediate A^(Hl) as compared 
to all clouds in the Hi complex. Cloud #20 has the lowest 
temperature and the smallest A^(Hl), and cloud #18 has an in- 
termediate temperature and the highest A^(Hl). The models, 
for hydrogen, are given in Table Note that t coo i > r co u for 
these representative clouds, indicating the validity of ioniza- 
tion equilibrium. 

6.3. Equilibrium Based On +5 and C 4 " 3 

For strict ionization equilibrium to hold, the ionization time 
scales must be smaller than the cooling time for all ions. 
The same arguments invoked for hydrogen above apply to the 
metals. Consider the +5 ion; the photoionization rate from 
the UVB at z = 0.67 is R ph = 4.3 x 10" 16 s" 1 , which yields 
T ph = 7 x 10 7 yr. 

For clouds #7+8, at T = 114,000 K, a r e C and a co u are 
3.7 x 10~ 12 and 3.8 x 10~ 16 cm 3 s" 1 , respectively. For logn„ = 
-4, we obtain r co u = r rec = 7 x 10 7 yr, where both scale in pro- 
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TABLE 4 

Cloud Models of the Hi Complex" 



Cld# 


jV(Hl) 


T 




fm 


a rec 


"coll 


''"cool 


Tcoll 






L 


M g 


L, 


Mj 


Note b 




[cm- 2 ] 


[K] 


[cnr 3 ] 




[cm 3 s _I ] 


[cm 3 s-'] 


[yr] 


[yr] 


[yr] 


[yr] 


[kpc] 


[M Q ] 


[kpc] 


[M Q ] 




/+o 


14. JU 


j.Ud 


— j.U 


— O.DZ 


— ij.Zl 


— o.zo 


o /in 


J. IK) 


1 n m 
1U.U1 


1 1 1 1 

11. /I 


A A 1 
4.43 


|C Q/| 

1 D.o4 


1 11 
L.I D 


i n "7i 
1U. 1 A 


T1 T 
J U 


7+8 


14.30 


5.06 


-4.0 


-5.69 


-13.21 


-8.28 


8.40 


4.70 


9.51 


9.78 


2.51 


11.05 


2.23 


10.22 


JU 


7+8 


14.30 


5.06 


-3.0 


-5.10 


-13.21 


-8.28 


7.40 


3.70 


9.01 


8.19 


0.91 


7.27 


1.73 


9.72 


EE 


1U 


1 3. /U 


A 1 1 
4. j 1 


— J.U 


a m 
— o.Uz 


1 O An 

— iz.OZ 


11 1A 
-1 1 Jo 




o. / J 


in ni 


lU.oo 


1 11 


1 1 11 
12. ZZ 


Z.JO 


Q A(\ 


J u 


10 


13.70 


4.31 


^t.o 


-5.02 


-12.62 


-11.36 


7.13 


1.16 


9.51 


8.88 


1.23 


7.22 


1.86 


9.10 


EE 


10 


13.70 


4.31 


-3.0 


^1.03 


-12.62 


-11.36 


6.14 


6.77 


9.01 


6.89 


-0.76 


2.25 


1.36 


8.60 


EE 


11 


13.43 


5.68 


-5.0 


-7.22 


-13.78 


-7.57 


10.86 


4.99 


10.01 


11.13 


4.17 


15.03 


3.04 


11.65 


JU 


11 


13.43 


5.68 


^1.0 


-6.48 


-13.78 


-7.57 


9.86 


3.99 


9.51 


9.39 


2.43 


10.81 


2.54 


11.15 


EE 


11 


13.43 


5.68 


-3.0 


-6.23 


-13.78 


-7.57 


8.86 


2.99 


9.01 


8.14 


1.17 


8.04 


2.04 


10.65 


EE 


15 


13.68 


4.95 


-5.0 


-6.53 


-13.12 


-8.50 


9.06 


5.92 


10.01 


11.05 


3.72 


13.69 


2.68 


10.56 


JU 


15 


13.68 


4.95 


^t.O 


-5.58 


-13.12 


-8.50 


8.06 


4.92 


9.51 


9.10 


1.77 


8.84 


2.18 


10.06 


EE 


15 


13.68 


4.95 


-3.0 


^1.88 


-13.12 


-8.50 


7.06 


3.93 


9.01 


7.40 


0.07 


4.73 


1.68 


9.56 


EE 


18 


14.26 


4.84 


-5.0 


-6.43 


-13.03 


-8.78 


9.06 


6.20 


10.01 


11.59 


4.20 


15.14 


2.62 


10.40 


JU 


18 


14.26 


4.84 


^1.0 


-5.46 


-13.03 


-8.78 


8.06 


5.20 


9.51 


9.62 


2.23 


10.23 


2.12 


9.90 


JU 


18 


14.26 


4.84 


-3.0 


-4.66 


-13.03 


-8.78 


7.06 


4.21 


9.01 


7.82 


0.43 


5.82 


1.62 


9.40 


EE 


20 


13.10 


4.41 


-5.0 


-6.09 


-12.69 


-10.64 


8.41 


8.05 


10.01 


10.30 


2.70 


10.63 


2.41 


9.75 


JU 


20 


13.10 


4.41 


^1.0 


-5.09 


-12.69 


-10.64 


7.41 


7.05 


9.51 


8.30 


0.70 


5.63 


1.91 


9.25 


EE 


20 


13.10 


4.41 


-3.0 


-4.11 


-12.69 


-10.64 


6.41 


6.07 


9.01 


6.32 


-1.28 


0.69 


1.41 


8.75 


EE 



a All quantities are base 10 logarithmic. The columns are: (1) cloud number; (2) hydrogen column density, logiV(Hl); (3) temperature, \ogT; 
(4) hydrogen density, log« H ; (5) hydrogen ionization fraction, log /hi; (6) recombination rate coefficient, logo rcc ; (7) collisional ionization rate 
coefficient, logo rec ; (8) cloud cooling timescale, logr coo i; (9) hydrogen collisional ionization time scale, logr co n; (10) cloud dynamical time scale, 
log 

r dym (11) cloud sound crossing time, logr sc ; (12) cloud absorption length scale, logL; (13) cloud gas mass, logM^; (14) cloud Jean's length, 
logL./, and (15) cloud Jean's gas mass, logMj. 

b JU denotes that the cloud is Jean's unstable or moderately Jean's unstable; EE denotes instability or slight instability toward expansion and/or 



evaporation. 

portion to n~ l . We thus have t co ii = r rec ~ 4r coo i independent 
of n„, indicating that the +5 ion is marginally in ionization 
equilibrium. Similarly, we find t cou ~ 3r coo i forC +3 . 

For cloud #10, at T = 20,000 K, a rec = 1.4 x 10" 11 cm 3 s" 1 
and a co n is negligible, indicating that collisional processes 
have vanishing importance in the ionization balance of +5 . 
For logn H = -5, we obtain r co n = T rec = 2 x 10 8 yr. For all den- 
sities, we have Tcooi/tcoII — 0.8. Similarly, for C +3 , we find 
Tbooi /t"co11 — 0.3. As with hydrogen, ionization equilibrium 
may be marginal for cloud #10 for O" 1 " 5 and C +3 . 

For cloud #11, at logT = 5.68, we have a rec = 1.2 x 10" 12 
and a co n = 6.9 x 10" 11 cm 3 s" 1 for O" 1 " 5 . When logn H = -4, 
and -5, we obtain r co n = 4 x 10 6-7 yr, where the longer times 
correspond to the lower density. We have r coo i / r co u ~ 2000 
for which O" 1 " 5 is predominantly in collisional equilibrium, i.e., 
Tboii/Tph — 20 if logn H = -4, and roughly equally balanced by 
photo and collisional equilibrium, i.e., t cou ~ r p h for logn H = 
-5. For C +3 , we obtain t coo \/t co \\ ~ 2 x 10 4 . Cloud #11 is 
clearly in ionization equilibrium. 

For clouds #15, #18, and #20, we find t C ooi/t C o11 — 2, 2, 
and 1, respectively for +5 . For C +3 , we find r coo i/r co n ~ 0.9, 
0.9, and 0.4, respectively. For these clouds, +5 is marginally 
in ionization equilibrium, whereas C +3 is marginally not in 
equilibrium. 

Note that the cooling times are 1-3 orders of magnitude 
shorter than the Hubble time, r H = \/Hq= 1.4 x 10 10 yr, ex- 
cept for cloud #11 for log n H < -5 (recall that the cooling time 
scales inversely with n H for fixed T). This would indicate that 
the cloud thermal conditions have evolved; however it is not 
possible to estimate the thermal histories of the clouds beyond 
speculating that they were hotter (and presumably less dense) 
at z > 0.67. For 4 < log T < 4.7 and for 5 < log T < 6, the 
cooling time increases toward higher temperatures, so we can 



infer that the thermal evolution of clouds in these tempera- 
ture was slower at epochs prior to z = 0.67 assuming the cloud 
densities have not strongly evolved. 

6.4. Dynamical Stability, Sizes, and Masses 

In view of the inferred thermal evolution of the absorbing 
clouds, it is interesting to examine their dynamical stability. 
This can be achieved by comparing the cloud dynamical and 
sound crossing time scales. Additionally, it would be of inter- 
est to estimate the cloud sizes and masses. 

To within a dimensiqnless factor of order unity, the dynam- 
ical time (see iLequeu x 2005, § 14.2.1), or free-fall time, of a 
cloud with total mass density p = p g /f g = (Auma/xufg)^, is 

Tdyn = 1 / y/Gp, where A H is the atomic weight of hydrogen, 
m a is the atomic mass unit, x s = 0.72 is the mass fraction of 
hydrogerQ and f g = ^It, /Vl m ~ 0. 16 is the fraction of the bary- 
onic mass in gas. Evaluating, we have Td yn ~ 3.2 x 10 7 n H ~ 1//2 
yr. 

The sound crossing time is r sc = L/c s , where L is the phys- 
ical length scale of the cloud and the sound speed is c 2 = 
jP g / p g = jkT I p,m a , where P g is the gas pressure, 7 = 5/3 for 
an ideal monatomic gas, and /i is the mean molecular weight. 
Our ionization models indicate that the clouds are highly ion- 
ized, so we assume p = (2x H + 3 / 4 x Hc + Vz-^z) -1 — 0.6, appro- 
priate for a fully ionized gas with logZ/Z Q ~ -2. We thus 
obtain, r sc ~ 6.4 x lO^LT -1 / 2 yr, when L is given in kpc. 

The "absorption" length scale of the cloud can be esti- 
mated as the line of sight path length required to give rise 
to the measured 7V(Hl) for the inferred n H from the ioniza- 
tion models, L = N(Hl)/{f m n H ). Assuming spherical geom- 

12 We have assumed a solar helium to hydrogen abundance ratio and so- 
lar metal abundance ratios for logZ/Zo = —2. A low metallicity primordial 
abundance pattern has Xn ~ 0.76, which yields a 5% difference in p s . 



12 



Churchill et al. 



etry, the gas mass of the cloud is crudely estimated from 
M g = p g L 3 ~ 3.4 x 10 7 n H i 3 M Q when L is given in kpc. A 
spherical geometry with unity volume filling factor is likely 
a very poor model, and one which will significantly overes- 
timate the cloud mass. If the absorber geometry is cylindri- 
cal of length L and radius R with aspect ratio j3 = 2R/L, then 
M cy i/M sp h oc 1 if the line of sight probes parallel to L. 

It is well established based on observational and theoretical 
grounds, t hat in general , Ly a forest clou ds cannot be pressure 
confined (Raucl fl998h . ISchayd J200l convincingly argues 
that clouds which develop due to the gravitational influence 
of underlying dark matter density perturbations persist in a 
near hydrodynamic equilibrium state local to the region giv- 
ing rise to the Lya absorption. Effectively, this is equivalent 
to stating that the dynamical time is equal to the sound cross- 
ing time, Tdy n ~ t sc , consequently implying that the character- 
istic length of the absorbin g region will be on the order of the 
local Jean's length, Lj (seelLequeux 20051 § 14.1.2). Equat- 
ing rdy n and t sc , the Jean's length is Ly ~ 5.0 x 10~ 3 (7 , /« H ) 1,/2 
kpc, from which the Jean's gas masQ can be estimated, 
My = p g L) ~ 4.3 T 3 / 2 ^" 1 / 2 M . 

For Mg > Mj the condition is Td yn < r sc ; the cloud is 
Jean's unstable to further gravitational contraction and will 
adjust on the dynamical time scale via fragmentation or due 
to shock processes. Conversely, for M g < Mj the condition 
is r sc < Tdy n , and the cloud will adjust in the s ound crossing 
time scale, eith er via evaporation or expansion (Schavel feOOll ; 
Lequeux 2005). The criterion assumes that the cloud is iso- 
lated, spherical, homogeneous, and exhibits no bulk motions. 

The dynamical and sound crossing times, deduced cloud 
sizes and cloud gas masses are listed in Table |4] for cloud 
#7+8, #10, #11, and the representative clouds #15, #18, and 
#20 for logn H = -3, -4, and -5. As described above, cloud 
#7+8 is well constrained to have log« H ~ -4, whereas cloud 
#10 and the remaining clouds with only limits on the metal 
line column densities are not well constrained. 

For log« H < -3, -4, and —5, we obtain T<j yn - 1 x 10 9 , 
3 x 10 9 , and 1 x 10 10 yr, respectively. Note that Td yn oc « H _1 , 
whereas r sc oc (/ HI n H ) _1 , and r coo i oc n H ~ l l 2 . The relative be- 
havior of these times scales is such that a transition from 
Tsc < T dyn to t sc > Td yn occurs as n H increases. This indi- 
cates that for larger n H , the clouds have a propensity to be 
in the regime of Jean's instability, whereas as for smaller 
n H , the clouds become unstable to evaporation and expan- 
sion. The last column of Table |4]notes whether a cloud model 
is inferred to be Jean's unstable (JU) or unstable to expan- 
sion and/or evaporation (EE). The behavior with n H suggests 
that between -3 < log« H < _ 5 there is a density at which the 
clouds would classify as being in hydrodynamic equilibrium. 
This is also reflected in the behavior of the Jean's lengths, Ly, 
which would be equivalent to the absorption length scales, L, 
when T dyn = 

Note that the absorption scale lengths, L, and the cloud gas 
masses, M g , become unphysically large as hydrogen density 
decreases. This would suggest that the clouds have logn H > 
-5, though we caution that the masses are probably overesti- 
mates due to the assumption of a spherical geometry. If more 
akin to cylindrical structures, viewed along the long axis, a 
factor of (3 = 0.1 reduces the masses by two orders of magni- 

13 The Jean's mass, Mj, usually applies to the total mass. If the gas frac- 
tion, f g = Qb/ Q,„ , in these clouds is near the cosmic mean, then gas mass and 
total mass are related by M tot = Mg/ f g . 



tude. However, the absorption length scale along the line of 
sight is independent of geometry. 

6.5. Hydrodynamic Equilibrium Conditions 

In view of the arguments given bv lSchavel (1200 lb that Ly a 
clouds should be in the regime of hydrostatic equilibrium, and 
in view of the unphysically large cloud scale lengths we de- 
duce in the regime of Jean's instability, it is reasonable to 
explore the inferred cloud conditions under the assumption 
of hydrodynamic equilibrium. This assumption will yield an 
equilibrium value for n H , from which the ionization conditions 
and metallicities can be constrained. 

Equating the dynamical time and the sound crossing time, 
and invoking the absorption scale length, L = A^(H I) /(f m n H ), 
we derive the condition of hydrodynamic equilibrium, 



n H x H f g jkT [/ H ,(« H , T)f = GA»pm\ [N(tti)f 



(2) 



where A^(Hl) and T are measured from the VP fit models to 
the Lya profiles, and f m (n H ,T) and n H are computed using 
ionization models of the clouds. 

Since the hydrogen ionization fraction is a function of n H , 
Eq. |2] is a transcendental equation and must therefore be 
solved numerically for the equilibrium hydrogen density for 
a cloud with temperature T. For each cloud, we interpolate 
a grid of ionization models with metallicity logZ m /Z Q = -2 
for —6.0 < logn„ < -2.0 in steps of Alogn„ = 0.2 and 2.6 < 
logT < 6.0 in steps of A log T =0.1. We set T equal to the 
measured cloud temperature (based upon the Doppler b pa- 
rameter from the VP fits), interpolate to obtain f m (n H ,T) as 
a function n H at that T, and iterate using Brent's method to 
locate the n H that satisfies Eq.|2]to a tolerance 1 x 10~ 10 . The 
method assumes that f m (n H ,T) is independent of metallicity 
(which we verified for log Z m /Z Q < -1). 

For this exercise, we treat clouds #7 and #8 separately (i.e., 
as a multiph ase s tructure as opposed to a single-phase cloud 
as done in § 15.11 ). We also omit cloud #9, which has a very 
uncertain Doppler b parameter in that its temperature is con- 
sistent with T = K. 

Once Eq. [2] is satisfied, the cloud absorption scale length, 
L, can be computed, from which the cloud gas mass can be 
estimated. Note that these quantities, while determined from 
the measured A^(Hl) and T as constrained by the ionization 
models, will be equivalent to the Jean's length and Jean's gas 
mass for the equilibrium n H . 

For each cloud, the uncertainty in the equilibrium n H and 
/ ra are obtained by accounting for the measured uncertainties 
in both Af(Hl) and T from the VP fits. The uncertainties in 
L account for the uncertainties in jV(Hl), n H , and f m . The 
uncertainties in M g account for the uncertainties in p g (due to 
n H ) and L. 

Once the equilibrium values are determined, the cloud 
metallicities and their uncertainties can be estimated from the 
equilibrium ionization cloud model. Denoting this metallicity 
as Z (in solar units), we have 



log Z = log 



N(X) 



iV(Hl) 



+ log 



Mn /in 



(3) 



eq 



where A^(Hl) is the measured neutral hydrogen column den- 
sity obtained from the VP fit to the data, N(X) is the measured 
or upper limit on the column density species X from the VP 
modeling, and n x is the number density of species X from the 
ionization model for the equilibrium f m and « H and the model 
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FIG. 8. — The deduced cloud physical conditions under the assumption of hydrodynamic equilibrium. — (a) The measured A'(Hl) from the VP fits. — (b) The 
measured temperature, based upon the VP Doppler b parameter. — (c) The equilibrium « H computed from Eq.|2] — (d) The equilibrium hydrogen ionization 
fraction, /hi, computed from Eq.|2] — (e) The metallicity, Z/Zq, computed from Eq.[5] the points are colored coded based upon the ion that provides the most 
stringent constraints; red for +5 , green for C +3 , and blue for C +2 . The dotted line is solar. — (/) The absorption scale length, L = N(Hl)/(f m n H ). Note that 
this is equivalent to the Jean's length, Lj, for the assumption of hydrodynamic equilibrium. — (g) The gas mass, M g . Note that M g = Mj, the Jean's gas mass, 
for the assumption of hydrodynamic equilibrium. — (K) The cooling time, t coo i; the dotted line is the Hubble time, t h = 1 /Hq. — (y The collisional time, t co h; 
the dotted line is r H . — (/) The dynamical time, fd yn ; the dotted line is t h . Note that the sound crossing time, t sc is equivalent to 7y yn for the assumption of 
hydrodynamic equilibrium. Open data points represent clouds that do not satisfy the criterion of hydrogen ionization equ ilibrium, where solid data points do 
satisfy the criterion. Dotted magenta boxes straddling clouds #7 and #8 represent the single-phase model presented in 8 15. li as illustrated in Figure [6] 



metallicity Z m . That is, we scale the equilibrium ionization 
model with Z m to obtain Z and an estimate of its uncertainty. 

In Figure [8] we present the deduced cloud properties un- 
der the assumption of hydrodynamic and ionization equilib- 
rium. Solid data points satisfy the condition of hydrogen ion- 
ization equilibrium, whereas open data points do not. The 
measured iV(Hl) and T are shown in panels [HJi and [8f>, re- 

fectively. The equilibrium n H and f m are presented in panels 
and[8jaf, respectively. Note that, as suggested by the analy- 
sis presented in Table|4] the clouds have -4.8 < log« H < -3.8 
and -6 < log/ UI < -4. 

The metallicities are plotted in Figure[8}?. Except for clouds 
#7, #8, and #10, the metallicities are upper limits. For the as- 
sumption of single ionization phase clouds, each metal line 
provides a unique limit. The most stringent limits are pre- 
sented, with the data point color coded by the ion that provides 
this best limit (redforO +5 , green for C +3 , and blue for C +2 ). In 
general, the upper limits are log Z/Z Q < —1.7. The measured 
values for clouds #7, #8, and #10 are logZ/Z Q = —1.651 
, respectively. 



-2.44!°;™, and- 



i 7O+0.20 



C+0.14 
-0.10' 



The absorption scale length (which is equal to the Jean's 
length under the assumption of hydrodynamic equilibrium), is 
shown in Figure[8}! for each cloud. The typical scale length is 
few hundred kpc, except for clouds #2, #3, #4, and #5, which 
have the lowest temperatures and thus relatively high densities 
and low ionization conditions. Recall, that the VP fits for 
clouds #1 through #6 should be viewed with caution. The 
gas masses, M g , of the clouds (which equal the Jean's gas 
masses under the assumption of hydrodynamic equilibrium) 
are presented in Figure |8^. Assuming spherical clouds, the 
gas masses are on the order of 10 8 -10 10 M Q , except for cloud 
#11, which has M g ~ 10 11 M©. These mass estimates should 
be viewed as upper limits by as much as one to two orders of 
magnitude. Cloud #1 1 is the highest temperature cloud with 
the highest ionization condition. The larger mass is a result of 
the large scale length and the fact that the hydrogen density is 
relatively high, logn H ~ —4.1. 

The cooling time, r coo i, is presented in panel [Sfi. Note 
that under the assumption of hydrodynamic equilibrium, the 
clouds have thermal stability in the order of 1-2 orders of 
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magnitude shorter than the Hubble time. The collisional time, 
Tcoii, is presented in panel|HI. Except for clouds #10 and #14, 
the clouds are in ionization equilibrium. However, clouds 
#2 through #6 also may not be in ionization equilibrium; we 
again remind the reader that the VP fits to these latter clouds 
are to be viewed with caution. The dynamical time (which is 
equal to the sound crossing time under the assumption of hy- 
drodynamic equilibrium), is presented in Figure [8]/' for each 
cloud. In all cases, Td yn and r sc are a factor of a few less than 
the Hubble time, which is indicated by the dashed line. 

The values obtained from this exercise are in remarkable 
agreement with t hose predicted from the simple scaling rela- 
tions proposed bv lSchavel (1200 ll) for the assumption of hydro- 
dynamic equilibrium. 

For comparison between the individual models o f clo uds #7 
and #8 and the combined cloud #7+8 modeled in § 15. U we in- 
dicate the results of the latter analysis as dashed boxes on Fig- 
ure [8] For cloud #7+8, the measured column densities were 
summed, from which the density and metallicity were simul- 
taneously constrained assuming a single ionization phase (see 
Figure 01. That analysis resulted in a slightly larger n a than 
the assumption of hydrodynamic equilibrium in the individual 
clouds, through the metallicities are consistent between anal- 
ysis methods. The difference in n a is likely due to the adding 
of the column densities. We also compare the individual cloud 
equilibrium and the combined cloud L, M g , t coo i, and Td yn t co u, 
shown as the dashed boxes on panels^- j. Note that cloud #7 
has Ten 3> T p h, where r p h = 10 4 yr, indicating that the hydro- 
gen in this cooler cloud is photoionized. On the other hand, 
cloud #8 has r co n ~ r p h for both hydrogen and +5 , and thus 
has a substantial collisional ionization contribution. 

Cloud #1 1 is of particular interest. This component is the 
hottest and most highly ionized, f m ~ 3 x 10~ 7 , cloud in the 
complex, and is among the highest density, logn,, ~ —4.1, of 
the clouds. We find r co n/r p h — 1 for hydrogen, and r co n/r p h < 
100 for +5 ; thus hydrogen is equally photo and collisionally 
ionized, whereas O" 1 " 5 , though not detected, is predominantly 
collisionally ionized. It is plausible that this cloud (#11) is 
shock heated gas, as further suggested by the fact that the de- 
duced cooling time (r coo i ~ 10 Gyr) is substantially longer 
than the deduced dynamical time (rd yn ~ 4 Gyr). Crudely 
adopting the dy namical time as a proxy for the compres- 
sion time (e.g., Birnboim & Dekel 2003; Dekel & Birnboim 
2006), and considering the temperature and ionization condi- 
tions, we find that cloud #1 1 is the only component in the Hi 
complex that is suggestive of shocked gas. 

6.6. Caveats 

The analysis we have presented has employed many sim- 
plifying assumptions. The VP fitting method philosophy as- 
sumes that the gas structure comprises several spatially dis- 
tinct isothermal clouds. In fact, it is very possible that the 
H I complex is a quasi-continuous nonuniform structure with 
temperature and density variations having a range of bulk mo- 
tions (and possibly at least one shock front, i.e., cloud #11). 
It is also plausible that such bulk motions can align in line of 
sight velocities creating caustics that emulate distinct clouds 
so that some VP components actually model a heterogeneous 
physical condition. 

Furthermore, the analysis invoking the dynamical time, 
sound crossing time, and Jean's mass and length is predicated 
on a homogeneous cloud. A cloud in thermodynamic equilib- 
rium cannot simultaneously be homogeneous and isothermal, 



as we have assumed here. 

In support of the assump tion of hydrodynamic equilibrium, 
we note that Schaye(2001) argues that a spherical cloud with 
an isothermal density profile, i.e., n H (r) oc r" 2 , has a well de- 
fined characteristic n H that is on the order of the maximum 
density probed by the line of sight. For the gas mass esti- 
mates, we have assumed spherical clouds, which is probably 
a very poor assumption. As such, the estimated gas masses 
should be considered upper limits. 

Finally, the ionization modeling assumes photoionization 
and collisional ionization equilibrium, which we have shown 
to be a valid condition for the majority, but not all of the 
clouds. In addition, the metallicity estimates are based upon 
the assumption of single phase ionization conditions. If some 
of the plausible concerns expressed above with regard to het- 
erogeneous physical conditions aligned in line of sight ve- 
locity hold, then multi-phase structure could be present that 
would affect the metallicity estimates. Overall, the assump- 
tion of ionization equilibrium in single phase gas is critical 
to all deduced quantities, especially the meta llicities and the 
thermal equilibrium values presented in § 16.51 

7. DISCUSSION 

With a velocity extent of 1600 km s" 1 , the Hi absorption 
complex at z = 0.672 in the quasar Q1317+277 is a most in- 
triguing gaseous structure. Absorption with velocity spreads 
on the order of 200 km s" 1 occur in approximately 10-15% 
of quasars (IWevmann et al.lfl99lt iGibson et al.ll2009l) . In al- 
most all cases, extreme absorption of this nature is produced 
by material ejected from the quasar itself (i.e., broad absorp- 
tion line quasars); however, the properties of the z = 0.672 
Hi complex studied here are not suggestive of absorption 
"associated" with or "intrinsic" to the quasar. For example, 
the Hi exhibits no evidence of part ial covering, which is an 
adopted signature of associated gas dBarlow & Sargent3 IT997l 
Gang uly et alJl 999). Furthermore, the kinematics of the met- 
als absorption lines are kinematical ly similar to the veloc- 
ity spreads observed i n galaxy halos (IChurchill & Vogtll2001l 
Churc hill et al.| [2003). Thus, the H I complex is likely iner- 
vening absorption. Early on, iBahcall & Salpeted (11965b sug- 
gested that the environments of galaxy clusters may give rise 
to extensive, intervening broad absorption line complexes, but 
few potential candidates have been identified. 

In this section, we summarize and further examine the na- 
ture and environment of the z = 0.672 Hi complex toward 
Ql 3 17+277, compare it to other similar Hi complexes, and 
discuss the possible origin of the Ql 3 17+277 Hi complex, 
such as hot-mode or cold-mode accretion, galactic winds, ac- 
creting filaments, intracluster gas, and/or the warm-hot phase 
of the intergalactic medium (WHIM). 

7.1. The Nature of the z = 0.672 Hi Complex 

The Q13 17+277 Hi complex at z = 0.672 is characterized 
by a velocity spread of 1600 km s" 1 and 21 Lya components 
with 12.9 < logTV(Hl) < 14.5, as determined by Voigt profile 
fitting. Under the assumption that the components (clouds) 
are near hydrodynamic equilibrium, the temperatures, hydro- 
gen number densities, and hydrogen ionization fractions pri- 
marily range between 10 4 < T < 10 5 K, -3.9 < log« H < -4.9, 
and -5.5 < \ogf m < -6.0. The deduced cloud sizes are on the 
order of 200 kpc, and the cloud baryonic gas masses range be- 
tween 10 9 -10 10 M©. Because the gas masses scale as n H f m , 
the large masses result from the low hydrogen number densi- 
ties and the high ionization conditions of the clouds. 
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The metallicities are measured only for clouds #7, #8, and 
#10 and are logZ/Z = —1.7, -2.5, and -1.8, respectively. 
Upper limits on the remaining clouds are log Z/Z < -2.4 to 
logZ/Z Q < —1.7. The limits on the cloud metallicities do not 
rule out enrichment at the level of the high red shift IGM (e.g., 
ICowie & Songailalll998l: iSimcoe et all 120041) . On the other 
hand, the low metallicities are 1-2 orders of magnitude below 
the log Z/ZfT) ~ -0.6 metallicities measured in z = 0.7 X-ray 
clusters dBalestra et al.ll2007l iMaughan et al.ll2008l) . 

Further insight is gained by examination of the kinematic- 
ionization substructure. Near the velocity center of the com- 
plex, at v ~ 75 km s -1 with respect to the galaxy Gl, is a hot 
T = 480,000 K, collisionally ionized component (cloud #11), 
which is very likely shocked gas. The narrow velocity re- 
gion just blueward of cloud #1 1 comprises four clouds, within 
which the only metal lines are detected. The general over- 
all absorption morphology of these four components (clouds 
#7-10), is that of a double profile suggesting two absorbing 
structures contiguous in velocity space. Cloud #10, separated 
by ~ 50 km s" 1 from cloud #11, has detected O VI absorption, 
which is deduced to arise in cool T = 20-30 K photoionized 
gas. Clouds #7 and #8, which have detected Ovi, Civ, and 
Cm absorption, give rise to a single profile, which is best 
modeled with a narrow core (cloud #7, cool photoionized gas 
with T ~ 10,000 K) and a broad component (cloud #8, hot 
collisionally ionized gas with T = 115,000 K). The velocity 
centroids are separated by less than half of a single COS spec- 
tral resolution element of ~ 17 km s" 1 . In addition, cloud #9 
appears to be a very narrow (T < 10,000 K) blue wing of 
cloud #10 offset by ~ 20 km s" 1 that overlaps with the red 
wing of cloud #8. These substructures may be suggestive of 
clouds moving through a hot, T > 10 6 K, medium in which a 
conductive interfaces aris es at the boundary b etween the cool, 
warm, and hot gas (e.g., Sembach et al. 2003|). 

Knowing the environment of the H I complex and relation- 
ship to galaxies would be instrumental for a broader interpre- 
tation. The proximate galaxy is Gl, which lies at D = 58 kpc 
from the quasar line of sight and has a redshift very near the 
mean of the Hi complex. The virial mass of galaxy Gl is es- 
timated to be with in a factor of two o f the virial mass of M87 
(M vil /M© ~ 10 14 . IStrader et al.ll2011l) . suggesting that galaxy 
Gl could be a central galaxy in a Virgo-like cluster. However, 
we find no clearly compelling evidence that Gl resides in a 
galaxy cluster or in a group with an X-ray emitting intraclus- 
ter medium. 

A search of the NASA Extragalactic Database (NED) and 
SIMBAD database yielded no reported X-ray measurements 
of Q1317+277. Within 4' of Q1317+277, there are no 
sources in the ROSAT all- sky survey bri ght source catalog 
dMickaelian et al.ll2006b IVoges et al.ll 1999t) . Of the five clos- 
est X-ray sources that are not identified either as a star or 
an AGN/quasar (which have known redshifts), one (1RXS 
J131954.5+253210)lies at 2° from Q1317+277 (50 Mpc pro- 
jected at z = 0.672). The bright R = 15.1 galaxy identified 
within 29 " of the X-rays would be 2 Mpc projected from this 
source at z = 0.672. If the X-ray source is associated with this 
galaxy, it is likely that the galaxy and X-ray source reside at a 
redshift much lower than the H I complex. 

In the ROSAT HR1 band, a minimum of 0.04 cnts s" 1 
is required for a source to be included in the ROSAT cat- 
alog. Applying th is upper limit, and invoking the relation- 
ship ((Mullis 2001) between the count rate and the total flux 
in the band (accounting for the aperture correction), we esti- 



mate L x < 2 x 10 39 erg s l , This is four orders of magnitude 
below the expected X-ray luminosity of 2 x 10 43 erg s -1 for a 
cluster with a central galaxy of virial mass of Gl, where we 
have employed th e bolometric X-ray l umino sity to virial mass 
scaling relation of Brvan & Norman (1998) and corrected for 
the X-ray band (see Mull isl 1200 ll) . The virial temperature of 
galaxy Gl is estimated to be on the order of 10 K, which 
yields a coronal tem perature o f kT c± 0.8 keV. According to 
the compilation of ICrain et al.l (1201 ll) . our upper limit on L x 
is not inconsistent with the observed X-ray luminosities of 
early-type galaxies with similar kT. 

We have successfully measured spectroscopic redshifts for 
only galaxies Gl and G2. As such, we cannot directly de- 
duce the presence of a cluster at z = 0.672, nor estimate the 
velocity dispersion of the galaxies that may reside at this red- 
shift. Based upon the photometric properties examined in 
the imaging date, it is difficult to definitively rule out or fa- 
vor the presence of a cluster at z = 0.67. However, the up- 
per limits on the X-ray flux within 50 Mpc projected from 
Ql 3 17+277 and the low metallicity of the H i complex, 1- 
2 dex below intracluster g as measurements (Balest raet al.l 
2007; Maughan et al. 2008), do not favor a large cluster nor a 
hot intracluster medium. 

Such considerations leave open the possibility that the Hi 
complex may be a phenomenon closely linked to a massive 
old elliptical galaxy that is not in an overdense environment. 
Taken together at face value, the data and the results of our 
analysis suggest a low metallicity structure, possibly a fila- 
ment or the remnants of a disrupted filament. There remains 
the question of the possible connection to the smaller galaxy 
G3, which unfortunately does not have a measured or esti- 
mated redshift. 

7.2. Review of Comparable Hi Complexes 

Given the dramatic velocity spread and kinematics of the 
Hi complex toward Q1317+277, and given its proximity to 
galaxy Gl (and possibly G3), it is of interest to investigate 
how rare/common are such absorbing complexes, what their 
observed relationships are with respect to galaxies, and what 
interpretations have been adopted in view of the role of gas in 
the evolution of individual and group galaxies. Such insights 
may help identify the Ql 3 17+277 Hi complex in a broader 
context. 

One example is the z — 2 Civ absorber complex toward 
the "Tololo Pair" (Tol 1037-27, z em = 2.18, and Tol 1038- 
27, z e .m = 2.33), whic h may be produced by intracluster gas 
( Jako bsen et all 119861) . The Civ doublets, later observed 
in two additional quasars in proximity on the sky, exhibit 
multiple discrete components with velocity widths ranging 
between 50-1000 km s" 1 and may extend some 18 Mpc 
dDinshaw & Impevll 19961) . 

Another possible intracluster absorption complex, at 
z = 0.695 toward the z em = 105 quasar PG 2302+029 
dJannuzi et al.ll 19961) . exhibits broad (Av = 3000 km s" 1 ) high 
ionization Civ, Nv, and Ovi doublets. In the FOS spectrum 
(velocity resolution Av = 230 kms" 1 ), the Lya absorption 
is segregated into multiple individual components each with 
Av < 250 km s" 1 distributed across the full velocity range of 
the metals. Near the central velocity, narrow Civ, Nv, and 
Ovi are present in one Lyq component. No low ionization 
species are present. Jan nuzi et al.l (119961) suggest three pos- 
sible interpretations: (1) material ejected from the quasar at 
extreme ejection velocity, (2) material associated with galax- 
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ies or the intracluster medium of a cluster or supercluster 
of galaxies, and (3) remnant material from supernovae in a 
galaxy. Unfortunately, their observations did not provide data 
capable of distinguishing between these scenarios; 

Tow ard the z em = 0.297 quasar H1821+643, Tripp et al. 
(120011) reported an Hi complex at z = 0.1212 comprising 
five Lya components distributed over a velocity interval of 
~ 700 km s" 1 in high resolution STIS and FUSE spectra. The 
logAf(Hl) column densities range from 12.7-13.8. Absorp- 
tion from O VI is present in a single broad wing of the central 
component, for which collisional ionization is favored with 
T = 10 5 3 " 5 6 K and -1.8 < [O/H] < -0.6. Seven galaxies in 
the velocity range of the absorption are present at impact pa- 
rameters ranging from 140-2400 kpc, with the 140 kpc galaxy 
aligne d in redshift with the Ovi absorption. Trip p et all 
(2001) favor the scenario in which the Hi complex is intra- 
group gas or an unvirialized filamentary structure. 

Using GHRS, STIS and FUSE spectra of the z em = 0.116BL 
Lac ob ject PKS 2155-304. IShull etalJ 0998) and lShull et af] 
(2003) analyzed an Hi complex with 14 Lya components 
centered at z = 0.056 spread over a velocity interval of 
2270 kms -1 . They estimate the Hi column densities have 
the range 14.5 < logN(Hl) < 15.0 with logZ/Z < -2.5 and 
cloud depths less than 400 kpc. Five H I emitting galaxies are 
found in the range 0.056 < z < 0.058 with impact parameters 
400-790 kpc. The two strongest Ly a blends have detectable 
Ovi and possible Ovill absorption measured in a Chandra 
spectrum. If the gas is the warm-hot ionized mediu m (WHIM) 
then t he density is constrained to logn H ~ -4. Shu ll et all 
(2003) favor a scenario in which the Ovi arises in "nearside" 
and "backside" shocked infall into the potential well of the 
galaxy group. 

In high resolution STIS an d FUSE spe c tra of the z em = 
0.370 quasar HS 0624+6907, lAracil et al.l (l2006al) report a 
cluster of 13 Lya lines[3 at z = 0.0635 with a velocity spread 
of 1000 kms" 1 . The Hi column densities range between 
12.6 < logAf(Hl) < 15.3. Only in the central component, with 
total logA^(Hl) = 15.4, are metal lines detected (Si III, Si IV, 
and Civ, but no Ovi) from which the gas is deduced to be 
photoionized with metallicity logZ/Z Q = -0.05, very near to 
solar enrichment, with log« H = -3.9. The gas temperatures are 
deduced to be T < 10 5 K. The estimated baryonic mass of this 
component is ~ 10 5 M Q with an absorption length scale of 3- 
5 kpc. They report 10 galaxies within 135-1370 kpc in the 
range 0.062 < z < 0.067, but this group is not consistent with 
elliptical-rich groups . On account of th e high metallicity and 
cool temperatures. lAracil et al.l d2006a) favor the interpreta- 
tion that this Hi complex is tidally stripped material from one 
of the nearby galaxies. 

The Hi complexes toward HI 82 1+643, PKS 2155-304, and 
HS 0624+6907 have both similarities and differences with the 
Hi complex toward Ql 3 17+277. However, the Q13 17+277 
Hi complex bares little resemblance to the metal-rich com- 
plexes observed toward PG 2303+029 and toward the Tol 
1037-27 and Tol 1038-27 pair. These latter two complexes 
may be examples of metal enriched intracluster gas. 

The broad Hi component in the complex toward 
HI 82 1+643 exhibits Ovi that is likely to be predominantly 
collisionally ionized with a relatively high metallicity. Simi- 
larly, cloud #8 in the Ql 3 17+277 Hi complex appears to be 

14 An "erratum" to Aracil et al. (2006a) was published (see Aracil et al. 
2006b). However, the deduced properties of the Hi absorbing complex are 
unaltered. 



a T > 10 5 K, collisionally ionized Ovi absorber, but accom- 
panied by C IV and C III absorption. On the other hand, the 
hottest, broad component in the Ql 3 17+277 Hi complex has 
no detected O VI and has upper limits on metallicity indicat- 
ing that it is metal poor in comparison. The HI 82 1+643 O VI 
absorber is at a substantially larger impact parameter to the 
nearest galaxy, which resides in a group that clearly has no 
massive elliptical galaxy, whereas the Q13 17+277 Hi com- 
plex quite is very close in projected to the massive elliptical 
galaxy Gl. 

The Hi column density for the low metallicity Ovi ab- 
sorber in the complex absorption toward PKS 2155-304 is 
1-2 orders of magnitude greater than the A^(Hl) of the com- 
ponents in the Ql 3 17+277 Hi complex. However, the clouds 
have similar n H . In both complexe s, the Ovi resides to the 
wings of largest Lya components. I Shull et all (|2003) inter- 
pret this as a shock interface, and this interpretation may ap- 
ply in the case of the Ql 3 17+277 Hi complex. However, as 
with the HI 82 1+643 Ovi absorber, the environment of the 
PKS 2155-304 Hi complex resides within a moderate group 
of galaxies having no massive elliptical galaxy. 

Of the three examples, the Hi complex toward HS 
0624+6907 has an Hi absorption profile morphology most 
similar to that of the Ql 3 17+277 Hi complex. The cool pho- 
toionized cloud with Civ, SilV, and Sim absorption com- 
pares to the cool photoionized Ovi absorbing cloud # 10 in 
the Ql 3 17+277 Hi complex, but cloud # 10 has no Si IV or 
C IV absorption. Though the clouds have similar n H , cloud #10 
has lower A^(Hl), higher ionization conditions, and a much 
lower metallicity as compared to the nearly solar metallicity 
cloud in the HS 0624+6907 complex. Because of the low ion- 
ization deduced for the latter cloud, the mass is ~ 4 orders 
of magnitude smaller than the mass deduced for cloud #10. 
Again, the galaxy environment of the HS 0624+6907 H I com- 
plex contains no massive elliptical galaxy. 

Overall, two unique features to the Ql 3 17+277 Hi complex 
are that it is at substantially higher redshift (z = 0.672) com- 
pared to the other reported Hi complexes (z = 0.056, 0.064 
and 0.121) and that it is in close projected proximity to a re- 
gion dominated by massive, metal-rich elliptical with an old 
stellar population. 

7.3. Interpreting the z = 0.672 Hi Complex 

Summarizing the gas properties of the Ql 3 17+277 Hi com- 
plex, we find (1) a hot, photo and collisionally ionized com- 
ponent that is consistent with shocked gas, (2) a cool compo- 
nent with photoionized O VI absorption, and (3) a cool compo- 
nent plausibly layered within a warm component that is both 
photo and collisionally ionized and exhibits Ovi, ClV and 
C III absorption, (4) several additional warm H I components 
spread over 1600 km s -1 in the rest-frame of z = 0.672, and 
(5) measurements of and limits on metal-enrichment between 
-2.5 <logZ/Z Q <-1.7. 

If the rest-frame velocity spread of the H I complex is due 
to the local Hubble flow, then the line-of-sight proper length 
of the structure would be D HF = Av/H^E(z), where Av = 1600 
kms" 1 and£ 2 (z) = fi m (l+z) 3 + Sl A . At z = 0.67, we estimate 
D HF ~ 15 Mpc. Based upon the ionization modeling, the de- 
duced physical sizes of the the H I complex are not consistent 
with a single structure of this extent. If the velocities of the 
Hi absorbers are due to Hubble flow, then the complex must 
comprise absorbers that are spatially segregated; there are five 
main absorption features apparent in the Ly a profile, which 
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would imply an average separation of 3 Mpc and that we have 
by chanced probed several isolated absorption systems. Con- 
sidering the similarities of the properties of the H I complexes 
toward H1821+643, PKS 2155-304, and HS 0624+6907, and 
that they are clearly connected with galaxies on scales of 0. 1- 
1 Mpc, we do not favor the interpretation that the Q 1 3 1 7+277 
Hi complex is multiple individual absorbers tracing a 15 Mpc 
Hubble flow. 

Is it reasonable that the individual peaks in the Hi column 
density are a result of having fragmented from a single struc - 
ture? If this structure were Jean's unstable, i.e., Ta yn -C r sc , 
then fragmentation and/or shock disruption into the observed 
components would be plausible. 

Crudely modeling this hypothetical single structure as hav- 
ing the Af(Hl) weighted mean temperature of the components, 
(logT) = 4.65 (prior to fragmentation or shock disruption, 
where the latter could subsequently heat the shocked com- 
ponents), and assuming that the length of its long axis is 
the sum of the Jean's lengths for each of the components, 
we find that Td yn ~ 72o T sc independent of the assumed mean 
hydrogen density, « H , where Td yn ~ 3 x 10 9 (n H /10~ 4 ) 1//2 and 
r sc ~ 6 x 10 I0 («h/10~ 4 ) 1//2 yr. The model is consistent with 
the plausibility that the multiple component structure of the 
Hi absorption complex could very well have resulted from the 
fragmentation of a quasi-coherent single gaseous structure. 

Our simple model yields a total Jean's length of Lj ~ 
2(« H /10" 4 ) 1 / 2 Mpc and a total Jean's mass of My ~ 2 x 
10 n (n H /10" 4 ) 1/2 Mpc, where the total mass is simply the 
sum of the Jean's masses for the individual components. The 
implied aspect ratio f3 = 2R/L for a cylindrical structure is 
(3 ~ (Mj/pgLj) ~ 0.1. If log« H = -4, then the Hi complex 
could be crudely envisioned as a 10 11 M Q mass cylindrical 
filament 2 Mpc long with a 200 kpc diameter. These values 
would scale as (n K ) l l 2 . 

Summarizing the luminous environment of the Hi com- 
plex, we find (1) a massive, red, high-metallicity, elliptical 
galaxy (Gl) with an old stellar population (~ 6 Gyr) residing 
at D = 58 kpc and aligned in velocity between the shocked Hi 
component and the metal-enriched Hi components, near the 
Hi profile velocity centroid, and (2) the upper limits on the 
X-ray luminosity consistent with the measured range for el- 
lipticals with kT ~ 0.8 keV and well below expected values 
for clusters dominated by massive ellipticals. 

In addition to the X-ray data not supporting the idea that 
galaxy Gl resides in a large cluster, the metallicity of the 
Hi absorbing gas is constrained to be 1-2 orders of mag- 
nitude below the average enrichment of clusters at z = 0.7 
dBalestra et all2007t Maugha n et al 112008b The observations 
of iBrinchmann & Ellisl (120001) and Ivan der Wei et al.1 (120051) 
indicate that old, massive, metal-rich late-type galaxies were 
present in the field as early as z = 2-3. The stellar population 
models for galaxy Gl suggest a formation epoch at z — 4. It 
may be that galaxy Gl is not a dominant elliptical galaxy of 
a large group, but may be a galaxy in the field that formed at 
high redshift. 

In addition to the aforementioned supporting observations, 
a massive, red, metal-rich elliptical far from a cluster envi- 
ronment i s theoretically plausibl e. The cosmological simu- 
lations of Gabor&Dave (2012), which treat accretion and 
feedback processes, yield a substantial fract ion of red, high 
mass g alaxies independent of overdensity. Johansson et al. 
(120 1 21) demonstrate that massive early-type galaxies are built 



in two phases. The first is an initial growth via in situ 
star formation fed by cold accret ion (z > 3) that is later 
quenched via vir i al shocking ("e.g.. iBirnboim & Deke 11 1 2003 ; 
i Keres et alJl2005blpekel & Birnboiml2006tlKeres et al] l2009 : 
Ivan de Voort et al.ll201 ll Ivan de Voort & Schaydl207lH . The 
second phase is accretion of stellar material (z ~ 1.5-0.5), 
called "dry minor mergers" due to the fact that their gas 
is heated and stripped by shock h eating in the hot virial 
halo of the massive galaxy (e.g., Khochfar&HjS [2009; 
i Hopkins et al.ll2010|). Th ese general results are also found by 
lNaabetal.1 (120071 120091) . 

With a mass of logM v i,/M Q > 13, galaxy Gl is well 
above the "critical mass" where the cooling time of 
the gas is much longer than the dynamical time, such 
that an a ccretion shock is established near the virial ra- 
dius (e.g.. IBirnboim & Dekelll2003HDekel & Birnboimll2006l: 
IKeres et alj 120091; Ivan de Voort & Schavd 1201 11) . The rate 
at which the cooling time increases with decreasing red- 
shift [r coo i oc l/p oc (1 + z)~ 3 ] is more rapid than the in- 
crease in the dynamical time with decreasing redshift, [rd yn oc 
1/y/poz (l+z)~ 3//2 ]. Thus, it is generally found in simu- 
lations that, at high redshifts, cold streams can often pen- 
etrate the hot halos and accrete onto the galaxies, whereas 
at intermediate to low redshifts, the longer cooling time re- 
sults in the heating and shocking of the cold fila ments, which 
then accrete into the halo but no t onto the galaxy (IKeres et al.1 
1 2005b iDekel & Bimboimll2006t iFaucher-Giguere et al.ll201 U 
Ivan de Voort et al.ll201 ll) . Galaxy Gl may have experienced 
the process in which it formed stars early during its initial 
accretion, perhaps even experiencing cold accretion from the 
filament that we are observing, and then at later times the ac- 
cretion penetrated no further than into the halo, resulting in a 
high mass, metal-rich, early-type galaxy with quenched star 
formation. 

Given such a scenario, the Hi complex could be interpreted 
as a filamentary structure with IGM chemical enrichment lev- 
els that is undergoing shock disruption near the viral radius of 
galaxy Gl. The hot cloud (#11) exhibits the signature of the 
shocked portion of the filament, and the multiphase absorbers 
(clouds #8-10) exhibit the signature of a conductive interface. 
Galaxy Gl is probed by the quasar sight line impact param- 
eter at D/R vlI ~ 0.1. Simulations indicate that "cold" fila- 
ments (i.e., those that have not been heated above log T = 5.5) 
penetrate no deeper than R ~ 0.5j? v i r at inter mediate redshifts 
(IKeres et"aTll2009t Ivan de Voort et al.1 120 1 ll) . suggesting that 
the sight line is intercepting the filament in the outer part of 
the virialized halo. 

A cartoon model of this interpretation is illustrated in Fig- 
ure [9] The cartoon is consistent with simulation results and 
is guided by inspection of Figure 7 (top panel, for a z = 2, 
Mvir = 10 12 M galaxy) from Ivan de Voort et all (1201 ll) and 
the "cold on ly" panel of Figure 6 for a z = 1 , M v ; r = 10 13 M Q 
galaxy from lKeres et al.l (120091) . 

Presumably, any components of the filament that have not 
yet been shocked will be shock heated and their cooling time 
will increase to be on the order of the Hubble time (as is 
deduced for cloud #11). Those components that may reside 
within the shock radius will likely not survive past R ~ 0.5/? v ; r 
and be assimilated into the hot halo. This would imply that, 
ultimately, the filament will accrete into the halo, be heated 
to the halo temperature, and never accrete onto the galaxy it- 
self. Galaxy Gl has likely not had new gas to fuel star forma- 
tion for several Gigayears and will likely not acquire new gas 
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Conductive Interface Region 
Stalled Shock Heated Gas 




FIG. 9. — A cartoon model illustrating the interpretation of the Hi absorbing complex as an accreting filament onto the galaxy Gl. The positive velocity material 
would require that it is infalling from the oberver's side of the galaxy, whereas the negative velocity material would be infalling opposite the observer's side of 
the galaxy. The plausible relative physical locations of the absorbing clouds, labeled by their VP component numbers (see Table|2j, are shown. The conductive 
interface (clouds #7-10) continue their infall after being stalled in a shock front, which is presently coincident with cloud #11 (see text). The model is adapted 
from simulation results (e.g., see the "cold only" panel of Figure 6 in Keres et al. 2009). 



via accretion as it evolves. Apart from the possibility of dry 
minor mergers dKhochfar & Silldl2009l iHopkins et al.ll20l0t 
iJohansson et alj|2012l) building the stellar mass of galaxy Gl, 
it is likely that the galaxy has experienced secular stellar evo- 
lution since z = 4. 

If galaxy G3 is at z — 0.67, one might ask if the filament is 
associated with this lower mass logM v j r /M ~ 11.9 galaxy. 
With the data in hand, it is nearly impossible to claim any re- 
lationship between galaxy G3, the Hi complex, and galaxy 
Gl. However, a virial shock is not predicted for lower mass 
galaxies, which have cooling times shorter than their dynami- 
cal times. As such, it would be expected that a filament accret- 
ing onto galaxy G3 would not show the signature of shocked 
gas with conductive interfaces as seen in the Hi complex. It 
could be speculated that galaxy G3 is embedded in the fila- 
ment, perhaps contributing to some in situ metal enrichment, 
and might eventually be a dry minor merger (2% mass ratio) 
with galaxy Gl. 

This scenario for the evolution of galaxy Gl is in stark 
contrast to the evolution of spiral galaxy G2, which is ob- 
served to have on-going star formation at z = 0.661, likely due 
to on-going multiphas e accretion that is feeding the galaxy 
dKacprzak et alj l2012h . To the limit of logJV(Hl) ~ 12.4, 
there is no indication in the COS spectrum of absorbing 
gas that bridges these two galaxies, which are separated by 
Av = 1960 km s" 1 (in the frame of the galaxies). 

The scenario of a shock heated filament for which the ac- 
cretion onto galaxy Gl is quenched is very different than the 
favored scenarios for the lower redshift H I compl exes toward 
H 1821 +643 PKS 2155-304, and HS 0624+6907 fShull et alJ 
[19981 120031: iTripp et al.ll200lt lAracil et al.ll2006al) . None of 
these complexes are in the vicinity of a massive elliptical 
galaxy, and as such it would be expected that they repre- 
sent physical scenarios other than the one we propose for 
the Ql 3 17+277 Hi complex. That is, Hi absorbing com- 
plexes, though filamentary in nature, likely trace various en- 
vironments around galaxies and their relationships with the 
intergalactic medium. 

8. CONCLUSIONS 



We have studied the gas properties of and the luminous en- 
vironment around a remarkable H I absorbing complex over 
the range z = 0.6686-0.6768 having a velocity spread of 
Av = 1600 km s" 1 towards the quasar Q1317+277. To con- 
strain the gas absorption properties, we have analyzed COS, 
STIS, and HIRES spectra of the quasar to examine the Ly- 
man series, Mgll, ClV, and Ovi absorption using Voigt pro- 
file fitting and photo+collisional ionization modeling. The 
galaxy properties in the quasar field have been measured using 
a WFPC2/F702W image, and multi-band APO/SPIcam and 
KPNO/IRIM ground-based images. Stellar population model- 
ing and halo abundance matching were employed to estimate 
the stellar ages, metallicities, mass, and virial mass, radii, and 
temperatures of the galaxies. 

Absorption Properties: The Hi complex is characterized by 
five contiguous absorbing regions, comprising 21 subcompo- 
nents, or clouds. Ionization models suggest that the kinemat- 
ics are not due to the local Hubble flow, but are consistent 
with a ~ 2 Mpc by ~ 200 kpc structure with a total gas mass 
on the order of 10 11 M©. We find a hot logT = 5.7 cloud, 
which we attribute to shock heated gas. This cloud is bor- 
dered by multiphase, cool, warm, and hot layers gas sugges- 
tive of a conductive interface in which metal line absorption 
(Cm, Civ, and Ovi) is detected. The low metallicity of the 
gas (-2.5 < logZ/Z Q < -1 .7) is more consistent with enrich- 
ment levels of the high redshift intergalactic medium (e.g., 
ICowie & Songailal [l998: Sim coe et al.l 2004). perhaps with 
low levels of in situ enrichment, than with th e metallicities 
observed in the intracluster medium at z = 0.7 fBalestra et al. 
l2007tlMaughan et al.ll2008h . 

Galaxy Properties: We find that the galaxy Gl (z = 0.6719 
at impact parameter D = 58 kpc) is consistent with a mas- 
sive (logM v i r /M Q ~ 13.7) elliptical galaxy that is metal-rich 
(Z ~ Zq), and formed at z = 4 (6 Gyr stellar population). The 
magnitudes and colors of the other galaxies in the field, limits 
on the X-ray luminosity, and the low metallicity of the Hi ab- 
sorbing gas are consistent with galaxy Gl being a field galaxy 
with an evolutionary scenario in which the cold accretion has 
been shut down by shock heating. The redshift of galaxy Gl 
(z = 0.6719) is bracketed by the redshifts of the collisionally 
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ionized gas (z = 0.6723, Av = +76 km s ') and the multiphase 
gas with detected metals consistent with a conductive inter- 
face (z = 0.6712-0.6717, Av = -124 to -30 km s" 1 ) in the Hi 
complex. 

Interpretation: Based upon our data and analysis, we fa- 
vor the scenario in which the Hi complex is a filamen- 
tary structure accreting into the halo of galaxy Gl that 
is experiencing virial shock heating and dynamical disrup- 
tion. Con sistent with predictions of both theoretical treat- 
ment (e.g.,[B irnboim jTDekell2003UDekel & Birnboiml2006t 
Birnboim et al. 2007j) and cosmological simulations (e.g., 
Keres et alJ2005Ll2009HKhochfar & Silkl2009HHopkins et alj 
2010t iF aucher-Gigu ere et al. 1201 It Ivan de Voort et al.l 1201 U 
van de Voort & Schavell201~ ; iJohansson et al.ll2012l) of mas- 
sive galaxies (i.e., logM v i r /M Q > 12), our observations and 
analysis indicate that the gas accreting into the halo of galaxy 
Gl at times after its formation epoch have not accreted onto 
the galaxy itself and that the Hi complex will also not accrete 
onto the galaxy. As such, the star formation of galaxy Gl has 
likely been quenched for Gigayears. 

Though it is difficult to definitively determine the nature of 
the Hi absorbing complex, the scenario we favor is highly 
plausible and consistent with simulations and theory. In fact, 
the data appear to provide convincing evidence that theory and 
simulations correctly predict by z < 1 , cold accretion via fila- 
ments in high mass galaxies is shock heated as it accreted into 
the virialized hot halos and that this gas then grows the halos, 
but does not fuel further star formation. Continued growth of 
the stellar mass would then necessarily occur via dry mergers 
of minor galaxies. 
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